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Abstract 

We prove the coarea formula for sufficiently smooth contact map- 
pings of Carnot manifolds. In particular, we investigate level surfaces 
of these mappings, and compare Riemannian and sub-Riemannian 
measures on them. Our main tool is the sharp asymptotic behavior 
of the Riemannian measure of the intersection of a tangent plane to a 
level surface and a sub-Riemannian ball. This calculation in particular 
implies that the sub-Riemannian measure of the set of characteristic 
points (i.e., the points at which the sub-Riemannian differential is 
degenerate) equals zero on almost every level set. 
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1 Introduction 



This article is devoted to a sub-Riemannian analog 

J J§ R {ip,x)dH v {x) = J dW{t) J dU v -\u) (1.1) 

M m <p-i( t ) 

of the well-known coarea formula (11.21) . Here p : M — > M is a smooth 
contact mapping of Carnot-Caratheodory spaces, and J~ R {(p, x) is a sub- 
Riemannian coarea factor determined by the values of a sub-Riemannian 
differential, or he- differential. In some sense, we may regard (ll.ip as a gen- 
eralization of the results of |Pa} IHej. Mgl , Mg2 and some other articles. 



It is well known that the coarea formula 

J J k (ip,x)dx = J dz J dH n ~ k (u), (1.2) 

where J^i^p^x) = ■>/ det (Dip(x)D(p*(x)), has many applications in analysis 
on Euclidean spaces [EH [FED EE21 ESI lOMSl ILTYal lot] . Here we assume 
that if e C\U, R k ) with U C W 1 for n > k. In particular, ([L2]) applies in the 
theory of exterior forms and currents, and in problems about minimal surfaces 
(see, for example, [FeFlj ) . Also, Stokes' formula easily follows from the coarea 
formula (see, for instance, [Vodlj ). The development of analysis on more 
general structures raises a natural question of extending the coarea formula 
to objects with a more general geometry than in Euclidean spaces, especially 
to metric spaces and sub-Riemannian manifolds. In 1999, L. Ambrosio and 
B. Kirchheim [AmKij proved an analog of the coarea formula for Lipschitz 
mappings defined on an "H"-rectifiable metric space with values in IR fc for 
n > k. In 2004, this formula was established for Lipschitz mappings defined 
on an % n -rectifiable metric space with values in an % fc -rectifiable metric space 
for n > k |Kmll IKm4] . Moreover, conditions were found on the image and 
preimage of a Lipschitz mapping defined on an % n -rectifiable metric space 
with values in an arbitrary metric space which are necessary and sufficient 
for the validity of the coarea formula. Independently of that result, the level 
sets of those mappings were investigated, and the metric analog of implicit 
function theorem was proved |Km21 IKm41 IKm6] . 

All results mentioned above were obtained on rectifiable metric spaces. 
Note that their metric structure is similar to that of Riemannian manifolds. 
However, there exist non-rectifiable metric spaces whose geometry is not com- 
parable to the Riemannian one. Carnot manifolds are of special interest (for 
results on their non-rectifiability, see, for instance, |AmKi] ) . Sub-Riemannian 
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geometry naturally arises in the theory of subelliptic equations, contact ge- 
ometry, optimal control theory, non-holonomic mechanics, neurobiology, and 
other areas (see, e. g., [NW1 l6ll ^ |Mgll IMTM^ 



IVod2[ IVod3[ ICS} IHlPaj ) . This theory has many applications. In addition, it 
has many well-known unsolved problems. 

One of them is the problem of the sub-Riemannian coarea formula, which 
is useful for developing a non-holonomic theory of currents, exterior forms, 
extremal surfaces (in sub-Riemannian and sub-Lorentzian geometries), and 
so forth. 

Heisenberg groups and Carnot groups are well-known particular cases of 
Carnot manifolds. In 1982, P. Pansu proved the coarea formula for real- 
valued functions defined on a Heisenberg group |Paj . Next, J. Heinonen 
|Hej extended this formula to smooth functions defined on a Carnot group. 
Another result concerning the analog of (11. 2p belongs to V. Magnani. In 



2000, he proved a coarea inequality for mappings of Carnot groups Mgl 
The equality was established for a mapping from a Heisenberg group to the 
Euclidean space M fc |Mg2|. The validity of coarea formula even for a model 



case of a mapping from a Carnot group to other Carnot group has remained 
an open question. 

The purpose of this paper is to prove the coarea formula for sufficiently 
smooth contact mappings (p : M — > M of Carnot manifolds. We emphasize 
that all results are new even in the particular case of a mapping of Carnot 
groups. 

As we mentioned above, for the first time a non-holonomic analogue of 
the coarea formula was discovered by P. Pansu |Paj . One of the basic ideas 
was to prove the sub-Riemannian coarea formula via the Riemannian one: 



V 



J§ R (m u) 

dU\z) I " , ' dH N - N (u) = I dW{z) I dU v -\u) 



1.3) 



Here N and N are the topological dimensions while v and v are the Hausdorff 
dimensions of the preimage and image respectively. It is clear that in sub- 
Riemannian spaces these dimensions differ. Other authors subsequently used 
the same idea. 

It follows easily from (11 .3p that the crucial point in this method is to 
understand the relation between the Riemannian and sub-Riemannian mea- 
sures on the Carnot manifolds themselves and on the level sets. Moreover, 
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an appropriate definition of the sub-Riemannian coarea factor is required. It 
is well known that the question of measure relation on the given spaces is 
trivial, whereas those of the geometry of level sets and the sub-Riemannian 
coarea factor are quite complicated. The main difficulties are related to the 
peculiarities of sub-Riemannian metrics. The inequivalence of Riemannian 
and sub-Riemannian metrics is illustrated for example by the fact that the 
Riemannian distance from the center of a radius r sub-Riemannian ball to its 
boundary varies from r to r M , M > 1, where the constant M depends on the 
structure of the space. Thus, an immediate question arises: how "sharply" 
does the tangent plane approximate a level set? (This question is important 
since the Riemannian tangency order o(r) is insufficient here: a level surface 
can "jump" out of a sub-Riemannian ball much earlier than necessary if this 
order of tangency is higher than o(r).) Also, it is a problem whether or not 
there exists a metric allowing us to describe the shape and geometry of the 
intersection of a sub-Riemannian ball and a level surface. Even if we answer 
these two questions, the most difficult question concerns the relation between 
the Hausdorff dimensions of the image and the intersection of a ball and a 
level set. 

In this article we solve problems stated above. Firstly, we divide the 
set of all points at which the classical differential is non- degenerate into 
a regular set and a characteristic set. Next, we define a sub-Riemannian 
quasimetric d 2 that enables us to calculate the measure of the intersection 
of a sub-Riemannian ball and a tangent plane (see below Remark 13.81 for 
details). A crucial idea in the construction of d<i is based on the fact that a 
ball in this quasimetric is asymptotically equal to direct product of Euclidean 
balls: 

B%(x,r) « B ni (x,r) x B n2 (x,r 2 ) x ... x B nM (x,r M ), M > 1, 

where N is the topological dimension of a Carnot manifold, and are the 
(topological) dimensions of the Euclidean balls B n \ % = 1, . . . , M. Thus, if 
a plane intersects a ball of this type then we can easily determine the shape 
of the intersection (contrary to the case of a "box" quasimetric with 

Box(x,r) ^Q ni (x,r) x Q n2 (x,r 2 ) x ... x Q nM (x,r M ), 

where we have Euclidean cubes Q ni , i = 1, . . . , M, instead of balls, since the 
cubes have sections of different shapes). Studying the sharpness of the ap- 
proximation to a level surface by its tangent plane, we introduce a "mixed" 
quasimetric possessing both Riemannian and sub-Riemannian properties. 
We prove that at regular points the tangent plane approximates the level sur- 
face sufficiently well, and this fact enables us to calculate the (Riemannian) 
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measure of this intersection, which depends on the Hausdorff dimensions of 
both the image and preimage. In other words, it is equivalent to r u ~ u (see 
Theorems 13 . 71 and 13 . 1 1 j) . These results yield analytic expression of the rela- 
tion between the Riemannian and sub-Riemannian measures at the regular 
points on level sets (see Theorem 13. 171) . 

The characteristic set case is a little more complicated since precisely 
near a characteristic point a surface can "jump" out of a sub-Riemannian 
ball. For this reason we cannot estimate the measure of the intersection 
of a ball and a surface via the measure of the intersection of a ball and the 
tangent plane at a characteristic point. Note also that in all papers mentioned 
above the preimage has a group structure, whose properties are essential for 
proving that the measure of the set of all characteristic points on each level 
surface equals zero. In the case of a mapping of two Carnot manifolds, both 
the image and preimage lack a group structure, and the approximation of 
a manifold by its local Carnot group is insufficient for the extension of the 
methods developed to our case. That is why we create an "intrinsic" method 
for studying the properties of the characteristic set. 

The result on the characteristic set is stated in Theorem 14.11 

In Section 4 we prove that the degenerate set of the classical differential 
does not affect either part of the coarea formula. 

Finally, in Section 5, we define the sub-Riemannian coarea factor (via 
the values of the foe-differential) and derive the sub-Riemannian coarea for- 
mula dUTj . 

For mappings if : H 1 — > R, formula (II. ip becomes 

\V HV {x)\- ^^dU u {x) = J dW{t) J dW~ 9 {u) 



since in this case we have 



N 



u A ■ 4 



Here Vh^(i) is the "horizontal part of the gradient of ip at x, and the 
coefficient appears due to the choice of a constant factor in the definition 
of the Hausdorff measure (to define W 1 we set this factor to u^, and if we 
consider 2 M instead of then this factor will be equal to 1); cf. |Paj . For 
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(p : HP ->■ R k , k < 2n, we deduce that 

f y/det(D H <p(x)D H <p(x)*) ■ U2n+1 ' U2n+2 ~ k dW{x) 

J ^2n+2 ' ^2n-k ' * 

M 

= J dU\t) J dU v -\u) 

M ¥> _1 (*) 

and 

jg*fo x) = y/det{D H <p(*)DH<p(x)*) • ^±L^±gzi , 

^2n+2 ' ^2n-k ' ^ 

Here Djj(p(x) is the "horizontal part" of the differential D<p(x), and if we 
replace all occurrences of Uu by 2 M then the coefficient following the Gram 



determinant of Djj(p(x) will be equal to 1; see also |Mg2|. Similarly, for 
M = G and M = R we obtain 

/ |V^(x)| U J ' Uu ~ l dW{x) = J dW{t) J dU v -\u) 



k=2 

and 



fc=2 

Once again, it is easy to see that if we take 2 M instead of a;^ then the coefficient 
following the length of the horizontal gradient of ip at x will be equal to 1; 
compare with [HeJ. 

Thus, for sufficiently smooth mappings we obtain all previous results as 
particular cases. 



2 Preliminaries 

In this section we introduce some necessary definitions and mention im- 
portant facts used to prove the main result. 

Definition 2.1 (cf. \Gr \ IKmVodt INSW] ). Fix a connected Riemannian C°°- 
manifold M of topological dimension N. The manifold M is called a Carnot- 
Caratheodory space if the tangent bundle TM has a filtration 

HM = HiM C . . . C HM C...C H M M = TM 
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by subbundles such that each point p G M has a neighborhood U C M 
equipped with a collection of (^'"-smooth vector fields Xi, . . . , X N , a G (0, 1], 
enjoying the following two properties. For each v G U, 

(1) HiM(v) = Hi{v) = span{Xi(u), . . . ,Xa m Hi{v)} is a subspace of T V M 
of a constant dimension dim H, i = 1, . . . , M; 

(2) we have 

[X i ,X j ](v)= c ijk {v)X k {v) (2.1) 

k: deg X j. <deg X^+deg Jf j 

where the degree degXk is defined as min{m | X k G H m }; 

Moreover, if the third condition holds then the Carnot-Caratheodory 
space is called the Carnot manifold: 

(3) the quotient mapping [-, • ] : Hi x Hj/Hj-i i— > Hj+i/Hj induced by 
the Lie bracket is an epimorphism for all 1 < j < M. 

The subbundle HM. is called horizontal. 

The number M is called the depth of the manifold M. 

Properties of Carnot-Caratheodory spaces and Carnot manifolds under 
assumptions of regularity mentioned in definition [2D] can be found in [KmVod[ 
IKm8l lKm9l IVodKm2l IKmlO] 

Example. A Carnot group is an example of a Carnot manifold. 
Definition 2.2. Consider the initial value problem 

N 

7 (t) = £^( 7 (t)), te [o,i], 
i=i 

.7(0) =x, 

where the vector fields X\, . . . , X^ are C 1 -smooth. Then, for the point y = 

( N \ 
7(1) we write y = expl ^ Vi^i ) W- 

,N s 

The mapping (yi, . . . , yjv) ^ expl y^Xj j (x) is called the exponential. 

\=i ' 

Definition 2.3. Consider 116M and (vi, . . . , vn) G Be(0, r), where Be(0, r) 
is a Euclidean ball in M . Define a mapping # u : Be(0,t) — > M as follows: 

N 

0„(vi, ... ,wat) = expl ^^Xj ) (it). 



r ) = expf y^f^ j 



It is known that 9 U is a C 1 -diffeomorphism if < r < r u for some r u > 0. 
The collection {vi}f =1 is called the normal coordinates or i/ie coordinates of 
the 1 st kind {with respect to u G M) of the point v = u (vi, ■ ■ ■ , vn)- 
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Theorem 2.4 ( |KmVodj ). Fix u G M. The coefficients 



Cijk{u) of ([271]) if degXi + degX,- = degX fc , 
otherwise 



define a graded nilpotent Lie algebra. 

We construct the Lie algebra q u of Theorem 12.41 as a graded nilpotent Lie 
algebra of vector fields {(X")'}^ on M such that the exponential mapping 

(xi, ...,Xiv) ^exp E W)' (°) is the identity p^lBLUj . In view of the 



8=1 7 

AT 



results of [FoSt j . the value of (XJ)'(O) is equal to a standard vector G 
where ij ^ ik if j ^ k, j = 1, . . . , N. We associate to each vector field of 
the resulting collection an index i so that {6 u )^{{X^)'){u) = Xi{u). By the 
construction, the vector fields {(X")'}^ satisfy 

mr,(xjy}= Yl ^wra' ( 2 - 2 ) 

deg X k =deg Xi+deg Xj 

everywhere on Mr. 

Notation 2.5. We use the following standard notation: for each X-dimen- 
sional multi-index \i = (//i, . . . , }1n), its homogeneous norm equals \fi\h = 

N 

Y^Hi degXi. 

8=1 

Definition 2.6. Refer as the graded nilpotent Carnot group G U M corre- 
sponding to the Lie algebra q u to the nilpotent tangent cone of M at u G M. 
We construct G U M in as a groupalgebra [Poj . that is, the exponential 
map is the identity: 

N 



exp(5]x i (X i ") , )(0) = (x 1 ,...,x JV ). 

»=i 



By the Baker-Campbell-Hausdorff formula, the group operation is defined so 
that the basis vector fields {X? )' on M N , i = l,...,N, are left-invariant |Po] : 
if 



N N 

8=1 i=l 



exp (J^ x t , y = exp (j^ Vi i?Q 



then 



N 

x ■ y = z = exp(j2 z^X?)'^. 



8=1 
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where 



z i = x i + y i , degXj = l, 

Zi = Xi + Vi + ^ K ,e 3 ( U ) ( X Wj ~ Vl X j ) > de § X i = 2 ' 
\ei+e 3 \ h =2, 

^ = ^+3/*+ KA^-y 13 ( 2 - 3 ) 

|ju+/8U=fc, 
At>0,^>0 

= Xi + yi + Gl c ^ l j (u)x tJ 'y p (xiy j - yiXj), degXi = k. 

\Li+e l +l3+e ] \ h =k, 
Kj 

Using the exponential mapping 9 U , we can push forward the vector fields 
(X?)' onto WcMas 

[(e u u(xry)}(e u (x)) = De u (x)((xn'(x)} 

and obtain the vector fields X? = Recall that X?(u) = X^u). 

Definition 2.7. Associated to the Lie algebra {Xf}f =l at u G M, is a local 
homogeneous group Q u Wl. Define it so that the mapping 9 U is a local group 
isomorphism between some neighborhoods of the identity elements of the 
groups G U M and G U M. 

The canonical Riemannian structure on £/"M is determined by the inner 
product at the identity element of Q U ~M coinciding with that on T U M. The 
canonical Riemannian structure on the nilpotent tangent cone G U M is defined 
so that the local group isomorphism 9 U is an isometry. 

Assumption 2.8. Henceforth we assume that the neighborhood U under 
consideration is such that U C Q u Ml for all u ElA. 

Definition 2.9. A curve 7 : [0, 1] — > M which is absolutely continuous in the 
Riemannian sense is called horizontal if ^(t) G ii/^M for almost all t G [0, 1] 
with respect to the Lebesgue measure on [0, 1]. 

A horizontal curve in G U ~M is defined similarly: here we require that 
7(t) G span{X 1 «(7(t)),...,^ d u im H 1 (7W)} = H" {t) M for almost all t G [0,1] 
with respect to the Lebesgue measure on [0,1]. 

Theorem 2.10 (see [21 EE] for smooth case and [KmVod[IKm9] for C 1,a -smooth 
case). Every two points o/M can be joined by a horizontal curve. 
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Definition 2.11. Given x, y G M, the Carnot-Caratheodory distance d cc (x, y) 
is defined as 

d cc (x,y) = inf{£( 7 ) : 7 : [0,1] M, j(t) G # tW M}, 

where the length £ of each (horizontal) curve 7 is calculated with respect to 
the Riemannian tensor on M. 

For x, y G 5"M the Carnot-Caratheodory distance d^ c (x, y) is defined as 

d u cc (x,y) = inf{r( 7 ) : 7 : [0, 1] ^ £"M, 7 (t) G ^ (t) M}, 

where the length £ u of each (horizontal) curve 7 is calculated with respect to 
the Riemannian tensor on Q U M. 

Denote the ball of radius r in d cc centered at x by B cc (x, r). Denote the 
ball of radius r in d u cc centered at x by B^ c (x, r). 

Assumption 2.12. Henceforth we assume that <p : M — > M is a mapping 
of Carnot manifolds M and M. We specify its smoothness below in Assump- 
tion E2SJ 

Notation 2.13. Hereinafter, we use the following notation. Namely, we 
denote: 

• the topological dimension of M (M) by iV (N); 

• the Hausdorff dimension of M (M) by v [y)\ 

• the depth of M (M) by M (M). 
In addition, we consider 

• horizontal subbundles H = #1 C TM and^H = H X C TM on M (M); 

• subspaces H C H 2 C . . . C H M = TM (H C H 2 C . . . C H M = TM) 
of dimensions n — n\ — dimi^! < dimH 2 < ... < dim Hm = N (h = 
h\ = dimifi < dimif2 < • • • < dimHj^ = N) respectively at every point 
x G M (x E M) (see Definition l2TT|) . 

We put rt\ = dimifi, h\ = dimifi, nk = dimiffc — dimiffc_i (n^ = 
dim H k - dim H k ^) z k = 2, . . . , M (M). 

Denote by d cc (d cc ) the Carnot-Caratheodory metric in M (M), and by 
d u cc (dj£) the Carnot-Caratheodory metric in ^"M (Q W M). 

Theorem 2.14 ( [Vod4] ). Suppose that E C M is an open set, and let <p : 
M — > M be a mapping with C 1 -smooth horizontal derivatives Xiip such that 
Xiip G H, i — 1, . . . , n. Then, it is hc-differentiable in points of E. Namely, 
for a point u G E, there exists a horizontal homomorphism L u : (Q U M, d" c ) — > 
(^^M, dtc ) of local Carnot groups such that it is continuous in u and 

d cc (tp(w), L u [w\) = o(d cc (u, w)) as E n Q U M 3 w — > u 

where o(-) is locally uniform. 
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Remark 2.15. Using the exponential mapping 6 U , we can consider L u both 
as a homomorphism of local Carnot groups and as a homomorphism of Lie 
algebras of these local Carnot groups. 

Notation 2.16. Henceforth we denote the /ic-differential L u of (p at u by 
the symbol D<p(u). 



Corollary 2.17 ( [Vod4] ). Let p : M ->■ M be a contact (i. e., D<p(H) C H) 
C 1 -mapping of Carnot manifolds (in the Riemannian sense). Then, it is 
continuously hc-differentiable everywhere on M. 

Property 2.18. Suppose that p G C^M^M). Then, the matrices of its 
differential Dcp(u) (in the bases {Xi(u)}f =l and {Xi{p{u))}f =l ) and of its 

/ic-differential D(p(u) (in the bases {X^{v)}f =l and {X i{p{v))}f =l , where v 
need not be equal to u) have the following structures: 



Dp(u) 



{ Vi(«) 










Vo(u) 



Vs(tx) 






/ Vii 



D<p(u) 



u 








V 



o 










Vo(u] 



VM o 








Vs?-i(«) 
















Vfir(«) 













* \ 

* 

* 

* 
* 



(2.4) 










0/ 



(2.5) 

where each block V, is of size (nj x n^), i = 1, . . . , M. Note that, the blocks 
Vi are the same in the matrices of Dtp(u) and Dip(u). 

Assumption 2.19. Throughout the article we assume that: 



N>N; 

If the matrix of D(p(x) has full rank at at least one point then we 
suppose that n > h. Otherwise, we suppose that 7ij > hi, i — 1, . . . , M 
(see Proposition I2.2ip . 
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Notation 2.20. Given yj <G C^M, M), denote by Z the set of points x G M 
with rank(D</?(x)) < N. 

Proposition 2.21. I. Given cp G C 1 (M, M), consider a point x G M \ Z . 

(a) 7/ify(aO(Vi) = then Dip(x){V{) = V { for alii = 1,...,M, and 
ia.nkD(p(x) = N. 

(b) IfD<f(x){Vi) C\ Vi then Y&nkD<p(x) < N. 

II. At the points x G M \ Z with rank D(p(x) = N we have 

(a) D<p(x)(H i ) = H i ,i = l,.„,M; _ 

(b) / D<f(x)(Hi/H i ^i) = Hi/H^i, i = 2, . . . , M. Here for every element 
y G Hi/Hi_i there exists Y G Hi such that 

y = {Y + T:Te = Y + H^, 

and we define for y the value T>ip(y) as the class 

Vip(y) = {Dip(Y) +Dip(T):Y G H t is fixed, T G H^} 

= Dip{Y)+Dip{H i _ 1 ). 

Proof. I. (a) Fix x G M \ Z such that D(p(x)(Vi) = V\ and denote D<p(x) 
by L. Verify that LV2 = V%. By taking into account property (3) of Definition 
12.11 and the property L[X, Y] = [LX,LY] of a group homomorphism L, 
we infer that L(V^) C V 2 . Verify that V 2 C Liy%). By property (4) of 
Definition 12.11 for each element Y G V 2 we have Y — \Yi, Y2}, where Y\, Y2 G 
V\. Since = LFj, where V, G Vi, z = 1, 2, it follows that 

Y = [Y h %\ = [LY X , LY 2 ] = L[Yi, Y 2 ] = LY, 

where Y G V 2 . Thus, D<p(x)(V 2 ) = %. 

Similar arguments show that Dip(x)(Vi) = Vi for all i — 3, . . . , M. Thus, 
D<p(x)[G x M\ = G v (a,)M, and rankLV(x) = AT. 

(b) Note that the image of every basis vector field under the /ic-differential 
Dip{x) is a vector field of the same degree or the zero vector field (i. e., 
the /ic-differential cannot decrease the degree of each basis vector field). If 
Dip{x)(V x ) £ V x then D^x)' 1 ^ \ [£ty>(x)(Vi)]) = 0, and D(p(x)[G x M] ^ 
G^wM. Consequently, rank Dip(x) < N. 

We prove claims (a) and (b) simultaneously. The properties of the Rie- 
mannian differential Dip(x) (see (12. 4p ) imply that D(p(x)(Hk) C H^ for 
k = 1, . . . , M. In particular, Dip(x)(Hi) C Hi by the definition of a contact 
mapping. Moreover, D(p(x)(Hi) = Hi if and only if D(p(x)(Vi) = V\. 
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In order to verify that Vip(x)(H 2 / Hi) C H 2 /H u take y E H 2 /H 1 . Con- 
sider y as a class of sums of a fixed element Y of H 2 , which may have nonzero 
coordinates (in the basis {Xj}^) only with indices greater than n\ and less 
than rt\ + n 2 + 1, and all elements of H\\ 

y = {Y + T :Y e H 2 is fixed, T e if J. 

Then, the structure of Dip(x) implies that the image 

= = D^(a;)y + D(p(x)(H) = Dip(x)Y + H e H 2 /H 1 

consists of the vectors with only the first hi + n 2 nonzero components. More- 
over, the components with indexes from hi + 1 to h 2 are the same for all 
vectors from y. (To verify that, it suffices to write the vectors Y and 
T in the basis {Xi}f =1 , and the matrix of Dtp in the bases {Xi}f =1 and 
{Xi\f =l1 and consider the product of this matrix and these vectors.) Thus, 
Vi P (x)(H 2 /H 1 )cH 2 /H 1 . 

Verify that Vip(x)(H 2 / H x ) D H 2 jH x . Take y e H 2 /H x and assume on 
the contrary that y £ T>(p(x)(H 2 / Hi). Consider a vector field F G J; then, 
Y = Y 2 + Yi with Y\ e H 1 . Without loss of generality we may assume that 
Yi = 0. Since rankD(p(x) = N, it follows that the images Dtp(x)(T x M.) and 
(D(p)-^(x)(T x M.) coincide, as in both cases we obtain the whole T^mM.. Here 
the linear operator (D(p)j^(x) acts on T x Wl (instead of G X M), and its matrix 

in the bases {X^x)}^ and {Xi(tp(x))}^L 1 equals that of Dip(x) (written in 

the bases {Xf(x)}f =1 and {Xp x \ip(x))}fL l ) . Denote by L the matrix of this 
operator. Then, there exists a vector Y with LY = Y. Without loss of 
generality we may assume that Y has at most n 2 nonzero components with 
indices from n\ + 1 to n\ + n 2 . Consequently, 

Y = LY = D<p(x)Y + [L- Dip(x)]Y , 

where Dtp(x)Y e H 2 and [L - Dtp(x)]Y e H 1 . Since Dtp(x)(H 1 ) = Hi, 
there exists Y x e D^(x) _1 ([L - Dcp(x)]Y ) n H v Put Y = Y + Y u and then 
Dip(x)Y = Y and Y G H 2 . Thus, Vip(x)(H 2 /Hi) = H 2 /Hi and, moreover, 
Dcp(x)(H 2 ) = H 2 . 

Similarly, we can show that T>ip(x)(Hi/ Hi-i) = Hi/H^ and D(p(x)(Hi) = 
Hi, % — 3, . . . , M. The proposition follows. □ 

Claim I implies 

Corollary 2.22. The existence of x e M with rank D<p(x) = N implies that 
rii > hi, i — 1, . . . , M . 
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Definition 2.23. The set 



X = {x G M \ Z : rank D<p(x) < N} 



is called the characteristic set. The points of x are called characteristic 
points. 

For t G M, denote the intersection </? -1 (£) H x by Xt- 

Remark 2.24. For M = R this definition coincides with the definition of 
the characteristic set given by P. Pansu |Paj and J. Heinonen [He]: in this 
case the condition rank Dtp(x) < 1 implies that all the horizontal derivatives 
Xi(p(x), % — 1, . . . , n, vanish at x, and vice versa. 

Remark 2.25. Proposition 12.211 explains the second part of Assumption 
12.191 if M 7^ x then we do not need to assume that rii > hi for i > 2. 
Assumption 2.26. If T-L v (x) — H N (x) = then we assume that ip G 



C^MjM) and X { e C 2 (M), i = 1,...,N (this condition is sufficient for 
establishing the /ic-differentiability of tp); otherwise, we assume that tp G 
C M+1 (M, M) and X { G C M+1 (M), i = 1, . . . , N. 

Definition 2.27. The set 



is called the regular set. If x G D then we say that x is a regular point. 
Lemma 2.28. I. For the set 



we have ( = (see Notation \2. 20\ for the description of Z). 

II. // there exists a family {X^, . . . , Xi~} of vector fields with the prop- 

e/rties rank ([Xi, <p] (x)) = N and deg X^ = v, then we have deg Xi, < M, 

3=1 



D = {x G M : rank D(p(x) = N} 




xeM\Z: 3{X h ,...,X if} } 



N 



(rank([A^](x))f =1 = N)=> (^degA i3 < u) } 




/ V\{u) * * * * 
V2(u) * * * 
V 3 (u) * * 



* \ 



* 



D<p(x) = 



* 












V o 







V 




* / 
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Choose N linearly independent columns with the minimal possible sum of 
the corresponding degrees (we consider degXj as the degree of column j, 



To this end, we must choose the maximal possible quantity of vectors 
from the blocks corresponding to the minimal degrees. In the first block, we 
can choose at most h\ linearly independent vectors. Next, take columns from 
rii + 1 to ni + n 2 and the corresponding "diagonal" block. In this block, we 
can choose at most h 2 linearly independent ^-dimensional elements. 

On assuming that there are more than h 2 linearly independent columns, 
we obtain a contradiction. Indeed, since the "diagonal" block is an [h 2 x n 2 )- 
matrix, there exists an elementary transformation reducing at least n 2 — h 2 
of its columns to zero. Apply this transformation to the columns of the 
matrix of Dip(x). Then, this block of size (hi + h 2 ) x (n\ + n 2 ) includes 
n\ + n 2 — h 2 columns of dimension h\ (more exactly, these column vectors 
belong to M™ 1 x O™ 2 ). Recall that the maximal number of linearly independent 
columns is hi, and we have already chosen them in the first "diagonal" block. 
Suppose that there are less than h\+h 2 linearly independent columns. Then, 
since rank D(p(x) = N, a "missing" column can be "replaced" by a column 
of a higher degree. Thus, we can obtain the minimal possible sum of degrees 
if we have only n\ — h\ + n 2 — h 2 linearly dependent among the first n\ + n 2 
columns. Therefore, the corresponding sum of degrees equals hi + 2h 2 . 

Applying further the same arguments to the degrees 3, . . . , M, we con- 
clude that the^minimal possible sum of degrees of linearly independent vector 
fields {X i:j {p}^ =1 is equal to v. Thus, claim I is proved. 

Claim II follows since^ firstly, we can obtain the sum equal to v only 
by considering the first M blocks, and, secondly if we have less than hk 
linearly independent vector fields on step k < M then the sum of degrees 
corresponding to the resulting collection is strictly greater than v. □ 

Theorem 2.29. I. The characteristic set x coincides with 



J 



1,...,N). 




N 




)} 



(2.6) 



II. The regular set D coincides with 



[xeM\Z: 3{X il ,...,X i .} 



N 



(rank([X^](x))f =1 = iV)^(^ 



deg Xi. = v 



)} 



(2.7) 
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Proof. I. (a) Denote the set in (12. 6p by A, and verify that A C x- Consider a 
point x G A. By (12. 6p . every N columns of the matrix of D<p{x) with indices 

N 

ii, . . . , ijy corresponding to a collection {X^, . . . , Xj^ } with degX^ = v 

are linearly dependent. Assume on the contrary that the rank of the matrix 
of D(p(x) equals N, and consequently, there exist N linearly independent 
columns in the matrix of Dip(x). The matrix of D(p(x) has a block structure, 
where the block k is an (n^ xnfc)-matrix. Thus, in each block, only columns 
can be linearly independent. Consequently, the sum of the degrees of the 
vector fields corresponding to these linearly independent columns equals v. 

The relation between the matrices of D(p(x) and D(p(x) implies that the 
corresponding columns of the matrix of D(p(x) are also linearly independent, 
and the sum of the degrees of the corresponding vector fields is equal to v. 
Thus, we arrive at a contradiction. 

The argument above implies that A C %. 

(b) Verify that x C A. Consider x G x- Since rank Dip(x) < N, it follows 
that every N columns are linearly dependent. Our goal is to show that if we 
take the columns of the matrix of Dip(x) with indices . . . such that 

N 

^2 deg X H - V 

then 

rank([X^](x))f =1 <iV. 

In view of claim I of Lemma I2.28[ it suffices to consider the columns with 
indices ii, . . . , ifj satisfying 

N 

^degJSQ. = v. 

3=1 

Take N columns Cj x , . . . , q_ of the matrix of Dip(x) corresponding to some 

N 

vectors X it , . . . , X ifj with ^ deg X^ = v. Since 

i=i 

rank({Q 1 ,...,c i ^}) < N, 

there exists a transformation r of the matrix of Dip(x) taking at least one 
of q. to a zero column, and preserving the block structure of the initial 
matrix. Denote this column by q . . Here we assume that jo is the minimal 
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number with this property. Put k = degX i:jQ — 1. Because of the block 
structure of the matrix Dtp(x), the numbers of nonzero entries of Cj iQ are 

at least dim ii^ + 1, and at most dimif^+i. Without loss of generality we 
may assume that ij = dim if^ + 1. Since we cannot transform all preceding 
columns into zero columns, it follows that j = dim if ^ + 1. 

Apply the same transformation r to the matrix of Dip(x), and consider 
the images of jo columns with indices ii, . . . , ij . Note that we may regard 
the image of column ij = dim + 1 as an element of IR dim Hk . Taking the 
structure of Dip(x) into account, we have j = dimH k + 1 vectors (columns) 
belonging to the space M dimHfc (because the components with indices greater 
than dim ii^ vanish). Thus, the rank of this collection equals dim if ^ < j . 
Consequently, the rank of {X^if, . . . , X if .<p} is strictly less than N. Since the 

collections of N vectors with the sum of degrees equal to v appear only in 
the matrix of Dip(x) (see Lemma \2. 281) . we see that if the sum of the degrees 
is equal to u, then the rank is strictly less than N. By taking Lemma 12.281 
and the fact that ia,nk(D(p(x)) = N into account, we see that \ C A. 

II. (a) Denote the set in (12.71) by B. Assume the contrary and take 
x G B with r&nk(Dip(x)) < N. Step I yields x G A, and thus, we obtain a 
contradiction since An B = 0. Consequently, BcD. 

(b) Take i6D. Since rank.D<£>(:r) = N, it follows that Y&nkDip(x) = N 
and x £ Z. Take an arbitrary collection of iV linearly independent columns 
of the matrix of Dip{x). Consequently, the corresponding columns of the 
matrix of Dip(x) are linearly independent as well, and x G B. Thus, DcB. 

The theorem follows. □ 

3 Properties of Level Sets 

In this section, we assume that x G M \ Z. 

First of all, we introduce a new metric, which is equivalent to the initial 
one and simplifies our computations. 

Definition 3.1. Let M be a Carnot manifold of topological dimension iV 

,N v 

and depth M, and put x = exp( %iXi J (g). Define the distance d 2 (x, g) as 



•i=i 
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follows: 

d 2 (x, g) = max j \ x j\ 2 ^j * 



1 

ni+ri2 i N 



(E E i*/)*" 1 "}- 



3=ni+l j=iV— nM+1 



A similar metric g?2 is introduced on the local Carnot group £/ u M. 
The set {y G M : ^(y, < r} is called the ball of radius r > centered 
at x and denoted by Box 2 (x,r). Similarly, Box%(x, r) stands for the ball in 



g?2 of radius r > centered at x. 

Remark 3.2. The preimage of Box 2 (x, r) in the metric d 2 under the mapping 
9 X equals 

Box 2 (0,r) = B% 1 (x,r) x B% 2 (x,r 2 ) x ... x B% M (x,r M ), 

where is a Euclidean ball of dimension rii, i — 1, . . . , M. 

Observe that Boxg (tt,r) = Box2(w,r) in the quasimetric d 2 since 



JV JV 

expl 



>(j2^r)(u) = exp(^^Xi)(u) (3.1) 

i=i i=i 

for all collections {x{\f =1 such that both parts of (13. II) make sense [Vod4j . 
The following proposition is useful for proving the main results. 

Proposition 3.3. Let M be a Carnot manifold of topological dimension N 
and depth M . Given a sufficiently small compact domain U (s M, there exist 
positive constants C > and r > depending on U , M , and N such that 
all points u and v of U satisfy 

[j Ba%(x,Z)CBaxZ(v,r + CO, < £, r < r . 

xGBoxj (v ,r) 

Proof. The proof follows the scheme of proof of the similar lemma for boxes 
in the metric of [VodKmlj IKmVodj . 

Put x = expf XiX^ ) (v), c?2(t>,a;) < r, and z = exp ( ZjX" J (a;), 
\=\ ' ' h=i ' 

d 2 (x,z) < £. Estimate the distance d^{v,z) applying group operation to 

points x and z. Namely, estimate the coefficients {Ci}£Li satisfying z = 
expfjt&xAiv). 
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Case of deg X { = 1. We have 

n n n n n 

5>) 2 = 5> + ^ = 5>) 2 + E(^) 2 + 2 ^ r + ^ 2desXi - 
i=i j=i «=i i=i j=i 

Case of degXj = 2. We have 

n+n 2 n+n 2 2 

J] (o) 2 < E + ^ + E K,eM)^ - z ^)) 

i=n+l i=n+l \ei+ej\ h =2, 

Kj 
n+n2 

< r 4 + £ 4 + 2r 2 £ 2 + E ( E ^,e,H(^- 



=n+l | ei + ej U=2, 



+ 2 E ((a* + iu){xiZj - ziXjij 

i=n+l \ei+ej\ h =2, 

< r 4 + £ 4 + Q(w)r 2 £ 2 + hiuyHj- + 2 

<(r + a 2 ( M )0 4 = (r + a 2 ( M )0 2desX % 

where bi(u) and q(-u) are linear combinations of the functions {F^ ej (u)}ij 
on assuming that \xiZj — z\Xj\ = 2r£ and \xi + z^\ = r + £ for all i,/, j. 
They are continuous with respect to u e [/. We can represent each sum 
r 4 + £ 4 + Cj(-u)r 2 £ 2 + bi(u)r^(r + £) 2 as (r + <ij(-u)£) 4 , where (ij(-u) depends 
on bi(u) and Put a 2 (u) = max rfj(-u) and assume without loss of 

i:degXi=2 

generality that (12(11) > 1. 

k 

Case of degXj = k > 2. Denote the sum n + n j by St- Then, as in 
the previous case of degXj = 2, we obtain 



E (^) 2 < E E 

lM+/8|h=fc, 
Ai>0,/3>0 

<(r + a fe ( M )0 2fe = (r + a fe ( M )0 2dcsXl . 



i=S fc _i+l i=5 fc _i+l |M+/?k=fe 

^>0,/3>0 



Here we use the property 



E |f;»|^-^< E 



fj,>0,P>0 m>0,/3>0 
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and define each function Ofe(u) in the similar way as a 2 {u) in the case deg Xi = 
2. We also assume without loss of generality that a k (u), p(u) > 1. Put 
a(u) = maxaj(w). The estimates above yield 



i 

<%(v,x) = max{(ElC/ X "' 

3=1 

n+n 2 i N ______ 

(E ior) 2dcsx " +i ,...,( E ioi 2 ) 2dcsXjv } 

j=n+l _=jV-n A/ + l 

< max{(r + aj(u)£) dc s x * } < r + a(tt)£. 

i 

Since all Oj(w) are continuous with respect to u, we may choose sufficiently 
large C < oo with a(w) < C for all w belonging to the given compact domain 
E/iM. The lemma follows. □ 

To prove the main theorems, we need a convenient quasimetric equivalent 
to the Riemannian metric. 

Definition 3.4. Given 

N 

" Yr) 



y = ex p (E ViX% 

8=1 



put p(_/,x) = max {|2/i|}. 

1=1, ...,jv 

Notation 3.5. In Theorems 13 . 71 and 13 . 1 1 \ we establish some local results for 
a fixed point x. From now on we use the auxiliary mapping ip = ip o 6 X . 

Notation 3.6. Put 



u (x) = min{ v : 3{X h , X } 



JV 

(rank([X^](x))f =1 = N)=> (E de g X - 

i=i 

It is clear that u \ x > v and z/ |d = 

Theorem 3.7. Fix x G <p _1 (t). Then, in a neighborhood in M . , the % N ~ N . 
measure o/7o[^> -1 (£)] flBox 2 (0,r) (see Remark \3.2ty is equal to 

Cr v - Va{x \l + o(l)) (3.2) 

where C is independent of r, and o(l) — )■ as r — >• 0. 
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Remark 3.8. We emphasize that the mapping if) — <p o 9 X acts on a neigh- 
borhood of the origin in M, N . Consequently, the tangent plane to the level set 
ifj~ l {t) lies in ~R N , and the intersection To[f/> -1 (t)] nBox2(0,r) is well-defined. 

Proof of Theorem \3. 7[ We split the proof into 6 steps. On step I we choose 
a suitable basis {wj}^ =1 for the tangent space T = T [^ _1 (t)] to the level set. 
On step II we define two projections of the basis vectors in T [?/> -1 (t)]. In 
particular, the first projection n assigns to each basis vector Wj some vector 
Pj of the same degree, and the second projection assigns a standard vector 
in {ex, . . . , e^v} to each basis vector Wj for T [^ _1 (t)]. On step III we 



show that rank(<7j , 0) 
vector Oi on if), i = 1, 



N 



= N (henceforth aiip stands for the action of the 
N), where 



{a u . . .,a~\ = {ei, ...,e N }\ 



71 



7Ti 



N-N 



N-N 



hi 



and on step IV we prove that the sum of the degrees of a iy % = 1, . . . , N, 
coincides with z/ (a;). Consequently, the sum of the degrees of itj, j = 
1, . . .,N — N, equals v — z/ (x). Further, on step V we deduce that the 
Lebesgue measure of Box 2 (0, r) flspanjp!, . . . ,p N _^} equals Cr u ~ u °( x \ where 
C is independent of r. Finally, on step VI we prove that the length of 
Mujj flBox2(0, r) equals 0(r k ^) for sufficiently small r > 0, and applying this 
result we show that 7r(Tn Box2(0,r)) coincides with the o(r)-neighborhood 
of S nBox 2 (0, r) in S, where S = span{p 1; . . . ,p N _j^} and o(r) is taken with 
respect to the metric d 2 . The theorem follows from the last result. 

Step I. Consider the normal coordinates at the point x. Recall that 
T = To[^ _1 (t)]. Choose an arbitrary basis in T, and write it as a matrix in 
which the basis vectors are written as rows. 

(i) By elementary row transformations, reduce this matrix to 



A 



* 
* 





* 

* 








\ 





V 









/ 



where the upper right triangle consists of zeroes, and the last nonzero entries 
of the rows appear in distinct columns. 

(ii) In the matrix A there is the following natural "grading" of the 
columns corresponding to the grading of the Lie algebra V of vector fields 



22 



on the local tangent cone. Split the columns of A into M blocks Bi, . . . , B M 
such that Bk consists of columns from dimi/ fc „ 1 + 1 to dimH k : 

1 • • • dim Hi dimifi+1 • • • dim_ff 2 dimif M -l+l • • ■ N 

( B. B 2 B M ). 

Next, there is also a "grading" of rows of A. There are M blocks 
Ai, . . . ,Am- Here A\ consists of the rows whose last nonzero element ap- 
pears in a column with index in [dim + 1, dim if;]: 

dimH2 . . . dimHjvj-i + 1 ... N \ 



... ... 

* * ... 



* 

* * * J 

Some of these blocks can be empty. 

(iii) For each k — 1, . . . , M, put Vk = Bk H Au 



/ Vi ... \ 

** V 2 ... ... 

, * * * ... 

; \ v. ; B-, \ v 2 * o 

* : : * Vm-i 



\ * * * * Bm-i \ Vm-i Vm ) 

Further, we transform the blocks Bk \ Vk, k = 1,...,M. Define the 
transformation of A by induction. For k = M, we have nothing to transform; 
thus, the base of induction holds. 

Suppose that we have transformed the blocks Bk \ Vk, k — I + 1, . . . , M, 
I < M — 1, and assume that Vi ^ (otherwise, we have nothing to trans- 
form). Replace the blocks Ai+i, ■ ■ ■ , Am by the projections of their row vec- 
tors onto (spanj^})- 1 H span{^4;, . . . , Am}- Roughly speaking, we remove 
the part collinear to Ai from the row vectors of Ai+i, . . . , Am- This projec- 
tion preserves the blocks Bk \ Vk, k — I + 1, . . . , M, I < M — 1, because of 
the "triangular" structure of A. 

Moreover, the rows of Bk \ Vk are orthogonal to Vk with respect to the 
classical inner product. 



A 2 



\ A M 



dim .Hi dimZfi+1 

* 
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Denote the resulting vectors by u>i, 



w 



N—N' 



( Vi 

* * V2 

* * * 



w 








\ 






v 9 J 



* 

V * 



V 

V 1 



M-l 



M-l 






V M 



( w 1 \ 



J 



\ W N-N I 



(3.3) 

Thus, we have constructed a "suitable" basis for T. 

Step II. For each j = 1, . . . , N — N, define the number k(j) as follows. 
Let be the row index of the last nonzero element of Wj in f)3.3p . and put 
Hj) = degX l{j) . 

Project each Wj onto VkU) by taking the last nonzero coordinates belong- 
ing to Vk(j), and denote this projection by pj = 7Tj-(wj). 

By the properties of W, this projection is orthogonal with respect to 



the standard Riemannian metric on 
(wj — 7Tj-(wj))l.Trj-(T)- Indeed, expand 



» N. 



for each vector Wj G T we have 



Wj 



Pj + Qj = ^r{wj) + qj 



and verify that qj is orthogonal to ir<r(T). Take u G itj(T). First, suppose 
that u is a "basis" vector in this image: u is one of the row-vectors of the 
block Vfc for some 1 < k < M. If pj G V; and / > k then qj±.u by choice of 
W. If I < k then obviously qj-Lu since qj and u have no nonzero coordinates 
with the same indices. Since j is an arbitrary number and u is an arbitrary 
basis vector, it follows that all qj are orthogonal to ttj-(T). Consequently, 
span{g,-, j = 1, . . . , N - N}±n T (T). 

Denote by iij the vector obtained from Wj by putting 



Ni(/)i 

for k ^ 




Put deg7Tj = k(j). By construction, the sum of the degrees of 7Ti, . 
equals the sum of degrees of the tangent vectors . . . , 

Step III. Consider the vectors 



• • • ' ^N-N 

)*( w n-n) ■ 



Wi, ■ ■ ■ ,o-jv} = i e i' • • • ,ejv} \ 



711 



7T 



7V-iV 



Kll' I rr 



N-N 



(3.4) 



of the standard basis in M. N . Since 9 X is a diffeomorphism, we have 

rank(([(6MU)WiM(x))L = N <=► ran^^]^))^ = iV. 
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Here the symbol (6 , x )*(0)((Tj) = D6 x (Q)(ai) stands for the action of the dif- 
ferential of 9 X at on the vector <jj, i = 1, . . . , N. 

Verify that the rank of ([o'iV'KO))^ equals N. Indeed, assume the con- 
trary Then ker Di/j(Q) n span{o"i, . . . , a^} ^ 0, and therefore 

dimT fl span{o"i, . . . , a^} > 1. 

Consider the coordinates of w G T fl span{<7i, . . . , cr^}. On the one hand, 
since w G T, the choice of the matrix W implies that we can expand w as 

N-N 
W = CLiWi. 
i=l 

Put i = max{i : a { ^ 0}. Since Wi = (w^ — 7Tj) + 7Tj, it follows that 

jV-iV 

w = gj(n;j-7rj) +a i Tr i 

i=l 

Recall that by the choice of {7Tj}^ JV , the only nonzero coordinate of Hi has 
index/ (io), and the nonzero coordinates of 

f y~] - ^i) + ^i] + ^0(^0 _ ^0)) 

have indices strictly less than l(io). On the other hand, since w G span{o"i, . . . , 
it follows from (13.41) that all coefficients of 7Ti, . . . ,7r N _ft in this expansion 
must vanish. Consequently, a io =0. Thus, we arrive at a contradiction. 

Step IV. Verify that the sum of degrees of basis vector fields {X^ }f = i 
such that Xi.(x) = (9 x )*(0)((jj), j = 1, . . . , N, equals vq(x). 

(i) Assume on the contrary that there exist standard vectors S\,...,S^ 
with 

N N 

M x ) = ^2 deg Si < deg Gi ( 3 - 5 ) 

i=l i=l 

and rank(5j , 0)ili = ^ (here deg<5j stands for the value deg[(9 x ) if (0)(Si)}, 
i = 1, . . . , N). Since {Si, . . . , 5$} 7^ {o"i, . . . , cr^}, there exists at least one 
= Sk for some j and k. Consequently, there exists a vector v satisfying 

v + rj^r G T. Moreover, the index of the last nonzero component of v is 
strictly less than /(j) (see the matrix W of (13.31) ). 



io-l 



'0 



7T 



so, 



T (2 j 7Tj . 
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(ii) Consider the chosen vector v. It is evident that v T (otherwise, 
Tij G T and <5 fc G T, which is impossible). Observe that rank(r l i/>)^! =1 < N, 
where = \{ v + ^~\) G 7", and r< = 5, for all i 7^ A;. Consequently, 

rank(A i '0)^ 1 = iV, where A^ = v and Aj = 5i for all i =fi k. Indeed, assume 
on the contrary that rank(Aj^)^ 1 < N. Then 

(6^)1, = 2(r^)f =1 - (A^)£i- (3.6) 
The assumption implies that f?/> = ^ai[$iV>] because rank(<5j'0)^ 1 = N. 

i=l 

Since 2t^ = (v + SyV' — 0, it follows from (13. 6p that the column of 

at 

{8iip)^ =1 is equal to — 1^/> = — X^ *^^]) an d, consequently, rank(5iV)i=i < 

i=i 

iV. Thus, we arrive at a contradiction. 

Since rank^^/Oili = N, assuming the contrary and applying similar 
arguments, we infer that there exists at least one coordinate (v)j of v such 

that rank(/v0)i=i = N, where )U& = ej and /ij = 5« for all i ^ k. 

(iii) Put 

Z = min{j : rank(//j^)^ 1 = iV", where = e^, /ij = 5i for z 7^ fe}. 

Note that I is strictly less than (the index of the only nonzero coordinate 
of Ti j). If dege; < deg7Tj then we obtain a contradiction with the assumption 

TV 

that Yl deg Si = vq(x). If dege; = deg^- then we consider the collection 

{Hi}f =l instead of {8i\f =11 and repeat the previous arguments. 

(iv) Observe that after finitely many iterations of the arguments de- 
scribed in substeps (ii)— (iii), we obtain a collection satisfying 

N N 

J^deg/ii < J^deg^ = u (x) 

i=l i=l 

(here deg/ij stands for the value deg[(d x )*(0) (jii)], i — 1, . . . , N). This con- 
tradicts to (13. 5p . 

Step V. (i) Denote by S the plane spanjpi, . . . , p N _^}, where the vectors 
Pi, . . . ,Vn-n are defined in Step II, and consider <SnBox 2 (0, r). Without loss 
of generality we may assume that Pi, ■ ■ ■ ,P N _ft are orthonormal (it suffices to 
carry out the procedure described in Step II within each V,, i = 1, . . . , M). 
Verify that this intersection equals a direct product of balls of radii r k , k = 
1, . . . , M, in some basis. 
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To this end, we consider a new orthonormal basis Oe 1: . . . , Oe N in M. N , 
where the orthogonal transformation O satisfies 0(span{Vj}) = span{Vi}, 
% = 1,...,M, and the image of each pj is the standard basis vector e^, 
j = 1, . . . , N — N. Then the (i 2 -bail in the initial basis coincides with the one 
in the new metric: 0(Box 2 (O, r)) = Box2(0,r). Thus, we have constructed a 
basis for S consisting of N — N standard basis vectors. 

Next, O(Box 2 (0,r) nS) = Box 2 (0,r) n 0(S) = Boxf (5) (0, r), where the 
latter is a ball in the metric d 2 \ ^ S y By the choice of O, it equals a direct 
product of Euclidean balls. 

C~)( *>\ 

(ii) Obviously, the (N— iV)-dimensional Lebesgue measure of Box 2 (0, r) 
equals Cr v ~ v °^ x \ where C depends only on x. Since O is an orthogonal trans- 
formation, so is the inverse mapping O^ 1 , and the (N — A r )-dimensional 
Lebesgue measure of Box 2 (0,r) n S — O~ 1 (Box 2 :)( - <s ^(0, r)) is the same; thus, 
it also equals Cr v ~ v ^ x \ 

Observe that at regular points the [N — A r )-dimensional Lebesgue measure 

M 

of S n Box 2 (0,r) equals f[ w nfc _ fifc r l/_z/ . 

k=i 

Step VI. On this step we show that there exists r > such that for 
r < r the set 

7r r (TnBox 2 (0,r)) 

coincides with the o(r )-neighborhood of S n Box 2 (0,r) in S, where o(r) is 
taken with respect to the metric d 2 - From this, we deduce that the % N ~ N - 
measure distortion of nj- equals 

Um ^(7r r (rnBox 2 (0,r))) = ^ ^(5nBox 2 (0,r)) 
■H N ~ N {mBox 2 {0,r)) »-^W JV - JV (rnBox 2 (0,r))' 

(i) Fix a vector Wj. Show that there exists r j > such that the length 
of Rwj n Box 2 (0,r) equals 0(r fc(j) ) for all r G (0,r 0j ). Indeed, fix r 0j > 
and a take point -u G M + u>j fl 9Box 2 (x, r j) and the corresponding vector u 
whose coordinates coincide with those of u. Here 



Wj = {z G R : z = awj, a G R+}. 



Then tt = Yl u k, where 

k=l 



u k = Uk(r j) e span{V fc }, 



and <i 2 (0,-u) = roj. Consequently, |ufe(roj)| < (roj) ; here | • | denotes the 
Euclidean length. 
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We verify that \u k (r)\ = 0{r k ^) for all r G (0,r 0j ), k = l,...,k(j). 
The vector u(r) = u k( r ) has the same coordinates as the point M + u>j PI 

k=l 

(9Box2(0, r). Assume on the contrary that for every roj and for every K < oo 
there exist r < roj and Ui(r) with 



|u,(r)| > Kr k( - j) . 

H 

|u(r)|- 



Fix such r , if and r < r , and put a r = rj\j- Then we can represent u(r) 



as 

fc(i) *(?) 
u(r) = a r u = ^a r u k = y\ fc (r). 
fc=i fc=i 

On the one hand, Kr k ^' < \ui(r)\ = \a r ui\ = \a r \\ui\. Therefore, 

r k(j) 

\a r \ > K- — -. 

N 

On the other hand, k = k(j), in view of its definition, satisfies \uk(j)(r)\ < 
LyH'j) f or a n r £ (o, r ]. Since Uk(j)(r) = a r ■ Uk(j), we have 



Lr k{j) > \u k{j) (r)\ = \ar\ • |u*(,-)| > ifr 



;,.j) kfcQ-)l 



Observe that depends only on Wj, and is independent of r > 0. Obvi- 

ously, the equality 

is violated for K > L. Thus we arrive at a contradiction. Consequently, 
\ Uk \ = 0(r fc W), r e (0, r ), for all Jfe = 1, ... , k(j). 

(ii) Consider Tfl Box 2 (0,r), and the projection tt-j- : T — > S. By the 
properties of the matrix W (see (13. 3p ). it is bijective. Verify that 

d 2 (0, d[n T (T n Box 2 (0, r))]) = r(l + o(l)), (3.7) 

where (9 stands for the boundary relative to the plane T. First, consider 
the basis elements w in T. Represent each of them as w = y + z, where 
7Tr {y + z) = y. The properties of ttj- yield deg z < deg y, where the degree of 

TV 

the vector u = u i e i is defined as 
i=i 

N 

degu = deg UjXj \ = max^jdegXj : Ui ^ 0}. 
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If w G T is not a basis vector, then we expand it in the basis {wi}^ =1 N as 



N-N N-N N-N 

w = } y a k w k = 2J afeZ/fc + 2j a k z k, 
k=i k=i k=i 

where 7rf(wk) — Vk, k — 1, . . . ,N — N. Next, it suffices to note that 

N-N N-N 
k=l k=l 

Thus, we can represent every weT&sw = y + z, where ir T (w) = y and 
deg y > deg z. 

Define a new quasimetric d 2E on Box 2 (0,r). For v,w G Box 2 (0,r) put 
d^ E (v , w) = 0^(0, it) — f ), where w — v denotes the Euclidean difference. This 
definition implies that Box2(0, r) coincides with the ball Box2_e(0, r) of radius 
r centered at in the metric d 2E . Consequently, for proving (13. 7p . it suffices 
to show that 

d° 2E (0, d[n r (T H Box 2 (0, r))]) = r(l + o(l)), 

where 9 stands for the boundary relative to the plane 7r<r(T). Take w = y + z 
with d 2 (0,w) = r and degy > deg 2;. Then 

d 2 (0 )2 / + ^) =d° 2E (0,y + z) =r, 

p(0,y + z) = 0(r dcgy ) (see Step V, substep (i)), and p(0,y) = 0(r dcgy ) since 
ttj- is a linear mapping. Then, we have p(y,y + z) = 0(r degy ) because yA-Z. 

I(z) 

Observe that d 2E (y,y + z) — d 2 (0,z). Represent z as z = z u where 

i=i 

Zi G Vj. Here degz/^) = 7(z) < degy because degz < degy. Since 
0(r degy )=p(y,y + z)~ max {|(^|}, 

l=l,...,I(Z) 
j=l,...,dim_ff i / J ff i _i 

it follows that |^| < 0(r degy ). 
This implies that 

1 1 
d 2 (0,z)= max < max {|zi| /(z) } 

i=l,...,I(z) i=l,...,I(z) 



< 0(r~*&>) < 0(r 1+ ^r) = r • o(l), 
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where o(l) is at most Cr M -^ for some < C < oo independent of the point 
of Box 2 (0,r). 

Proposition^. 3limplies that d% E (v,w) < d 2E ( f, u)-\-cd^ E {u, w) for u,v,w G 
Box 2 (0,r) and some c > (it suffices to put F l p(x) = in Proposition 13. 3p . 
Consequently, 

d 2 (0, y) = d° B (0, y) < d° 2E (0, z + y) + cd° 2E (y, z + y) = r(l + o(l)), 
and (13. 7p follows. 

Thus, the measure of Tfl Box 2 (x, r) is equivalent to Cr v ~ uo ^ x ' as r — > 0. 
The theorem follows. □ 

Definition 3.9. Let £ : M — )■ M be a mapping of two Carnot manifolds. Fix 

x G M. The d 2 - distortion of £ at y ifit/i respect to x equals 

^2(3/, a;) 

and the p-distortion equals 

p(y,x) 

To simplify notation, we denote the Gram determinant y^det{AA*) for 
a matrix A by T>(A). We denote the Gram determinant a/ det(5*I?) for a 
matrix 5 by T>(B). 

Proposition 3.10. The matrix of the differential of the mapping 6 X , x G M ; 
at equals the identity matrix. 

Henceforth we denote the Riemann tensor at y by g(y). 

Theorem 3.11. Suppose that x G y3 _1 (t) is a regular point. Then 

(I) In a neighborhood ofO = 9~ l (x) there exists a mapping from To[^ _1 (t)]fl 
Box 2 (0, r(l + o(l))) to ip ~ 1 {t) flBox 2 (0, r) such that both d 2 - and p- distortions 
with respect to are equal to 1 + o(l) at every y G To["0 -1 (t)] D Box 2 (0, r(l + 
o(l))), where o(l) zs uniform in x = 6 X (0) GWdM ana 1 in y & T [ifj~ l (t)] (1 
Box 2 (0,r(l + o(l))J. 

(II) TTie H N ~ N -measure o/</? _1 (t) fl Box 2 (x, r) equals 

V(g\ ke r Dv{x) ) ■ fl ^n fc -n, • ^^ ^(l + o(l)) 
fc=l £>(£V(x)) 

with o(l) —7-0 as r — > 0, where o(l) — >■ uniformly in x G W <e M. 
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In the proof of Theorem 13 .111 we use the following notation introduced in 
Theorem 13 .71 We denote the mapping tp o 9 X by ip, and the tangent space 
T ['?/>~ 1 (£)] by T. We also use the mapping ttj- defined in Step III and the 
image S = ttj-(T). 

Proof of Theorem \3.11[ Notice that without loss of generality we may as- 
sume that Dip(z) is strictly separated from on (ker DiplO)) 1 - for z G Box 2 (0, r), 
i. e., Di[)(z)(v) > a > for all z G Box2(0,r) and v G (ker Dip^O)) 1 - with 
\v\ = 1 (here | • | denotes the Euclidean norm). 

Step I. Verify that C in (13.21) at a regular point equals up to a 

Riemannian factor. In particular, we show that ^^M x ^ equals the measure 

1 T>{D<p{x)) 1 



distortion under the mapping ttj- defined in Step II of Theorem 13.71 

To this end, consider the projection of the normal space M to T onto 
H^(0) constructed by the choice of vectors written as rows in the matrix of 
Dip(0) (recall that the matrix of Dtp(0) equals the matrix of ([(6'~ 1 )*X 1 ]'?/>, . . . , [(6 , ~ 1 )*Xjv]V')> 
and in the bases {e i }^ 1 and {e i }^ 1 the matrices of Dtp(0) and Dip(0) have 
the structure similar to (I2.4p and (12.51) . respectively). Denote this projection 
by ii^f and its image by £ = nj^(Af). Verify that S is orthogonal to C. Write 
a basis for Af as the row vectors of the matrix Dip(0) in the basis {ei\f =1 . 

(i) Without loss of generality we may assume that ttj^ is an orthogonal 
projection (see Theorem 13. 7p . Indeed, the matrix of Dip(0) admits "gradings" 
Bi, ... , Bj-f of the columns and Ai, . . . , of the rows similar to those of 
the matrix A of Step II substep (ii) of Theorem 13.71 Namely, we say that 
a column vector belongs to the block Bk if its index is at least dimiffc_i + 1 
(here we assume that dimifo — 0) and at most dimiffc, k = 1, . . . , M, and 
we say that a row vector belongs to the block Ai if the index of its first 
nonzero element is at least dimiffc_! + 1 (here we assume that dimiJ = 0) 
and at most dim Hk, k = 1, . . . , M^Put Vk = Bk H Ak^k = 1, . . . , M. 

Next, we transform the blocks Bj. \ k — 1, . . . , M. Define the required 
transformation of the matrix of Dip(0) by induction. For k = 1, we have 
nothing to transform; thus, the base of induction holds. Suppose that we have 
transformed the blocks Bk \ 14, k — 1, . . . , I, I > 1. Assume that Vi + \ ^ 
(otherwise, the transformation is trivial). Replace the blocks Ai, . . . ,Ai by 
the projections of their row vectors onto the space 

(span^+i})- 1 n span{^.i, . . . , Ai +1 }. 

This projection leaves the blocks Bk \ 14, k = 1, . . . I > 1, unchanged 
because of the "triangular" structure of the matrix of Dip(0). 
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Moreover, the rows of the block Bk \ Vk are orthogonal to those of Vk, 
k = 1 M. It results 
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By the construction, the blocks Vi, . . . , Vj^ constitute the matrix of Dip(0) 
(see (I2.5p ). and dim£ = N at regular points. 

(ii) It is easy to see that C is orthogonal to S. 

Indeed, take two vectors v G C and w G S. First, suppose that v and w 
are the images of arbitrary basis vectors in Af and T respectively. Then 
v G Vi and iy G Vj for some i and j. Hi ^ j then obviously V-Lw (since they 
have no corresponding nonzero components). 

Suppose that i — j and consider the preimages of v and w. By construc- 
tion, they are orthogonal. Note that the preimage of v belongs to A/", and 
we can write it as (0, . . . , 0, v, *), where the nonzero part begins with some 
element v of v. The preimage of w belongs to T, and it can be written as 
(*, w, 0, . . . , 0), where the last component of the nonzero part is some element 
w of w. Consequently, 

= (Q,...,Q,v,*)-(*,w,0,...,0) T = (v,w), 

and V-Lw. 

Since v and w are the images of arbitrary basis vectors in T and M, we 
have S-LC 

(iii) Thus, S is orthogonal to C. Moreover, the projection of T onto 
S is also orthogonal (see Theorem 13. 7ft : (v — n-j-(v))-Ln T (v) for all v G T. 
Let us verify that the measure distortions under the projections ttj- and n^f 
coincide. 

Step II. (i) Show that 

dimT-dim(T H S) = dim A/" - dim(A/" n £). 

Put / = dim(T PI S). On the one hand, since A/"_LT and £_LiS, it follows 
that span{A/",£}±(T nS), and so span{A/",£} C (T fl^) 1 . On the other 
hand, if i>_L(T D S) then t> G span{A^, £}. Indeed, assume on the contrary 
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that v span{A/", £}. Then v = vmc + v (M £)■*-■> where vj^c £ span{A/", £} and 
V {NC) 1 - e (span{A/", -C})- 1 with v^ C )± ^ 0. Consequently, since v^ C )±-LAf 
and W(_^/-£)±_L£, we obtain vos£\± G (Tfl <S) and arrive at a contradiction. 

Thus, span{A/", £} = (T fl <S) ± , and dim span {A/", £} = N — I. Since 
dim A/" = dim£, we get 

dim(A/ r n £) = 2 dim Af - dim span{A/", C} = 2N - N + I. 

Consequently, 

dimT-dim(TnS) = N-N-l = dim A/" - dim(AT n C). (3.8) 

(ii) Recall that £ = kerZ>0(O) at regular points (see Step I, substep (ii)). 
Verify that dim(T Pi C) — dim(A/" fl S) = 0. Indeed, assume on the contrary 
that there exists h G TC\jC, h ^ 0. Then, h-LAf and hJS. Since (Tn/2)- 1 = 
span{A/", <S}, which we can justify similarly to the substep (i) above, it follows 
that dimspan{A^, S} < N — 1. Consequently, from dimspan{A/"} = N and 
dimspan{5} = N — N, we infer that dim(A/" fl «S) > 1 and there exists 
h 1 - G (AfdS); thus, Z>0(O) vanishes on some vector in Af. This implies that 
one of N vectors in Af is orthogonal to C Consequently, its image under 
iij^f vanishes since every vector in Af can be uniquely represented as a sum 
of a vector in C and a vector in C L . Thus, ir^ has a nonzero kernel, and 
dim(7TA/-(A0) < N- This implies that dim [7^ (AO] = dim[(ker i3^(0))- L ] < N, 
which contradicts the regularity of x. 

Hence, dim(T n C) = and dim(M n S) = 0. 

(iii) Put g = dimT — dim(T fl S). Consider the maximal subset S' of 
S satisfying S'±(T H 5). Namely, 5' = (Tn «S) ± n 5. Then dim 5' = g. 
Put T' = 7r 7 :1 ( l S / ). It is easy to see that T'l(Tn 5). Indeed, take arbitrary 
vectors i> G T fl <S, i>' G S' , and a vector u>, such that t> ' + w = Tiy 1 ^'). In 
particular, v' + w G T. Note that u>_L(T H 5) because 7rr is an orthogonal 
projection. Then T 9 itZ^(y') = (t>' + ty)_Lv. Since v and t>' are arbitrary, it 
follows that T'-L(Tn«S). The non-degeneracy of 717- and tt^- 1 yields dimT' = 

Thus, we may regard the following mapping as the "inverse" projection 
-Kj- 1 . To each v G S it assigns the vector iu G span{v, £} fl T such that 
w = (!) + £)nT. It is well defined since £_LiS and dim(T fl C) — 0. 

(iv) Since £ fl A/" is orthogonal to both <S and T, we have 

(v + C)nr = (v + S)nr 

for v & S, where £ C (£ fl A/") -1 fl £ is the minimal subset enjoying this 
property. Let us prove that £ = (C fl Af) 1 - fl £. It is easy to see that if 
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w = (v+£)(~)T = v+u with u E £ then u = w — v G span{5, T} = (CnAf) 1 . 
Consequently, «G(£n Af) L fl £, and dim£ < N — N — I. 

Verify that dim£ > N — N — I. Assume on the contrary that £ 7^ 
(£rW)- L n£, and consider w G £:- L n((£nA/')- L n£). Since v G span{T,S} = 
(£ fl A/") -1 , it follows that vo = vj + Vq , where vJeT and Vq G S. Take 
v = — G <S and consider v q = v + v G T . By the assumption on £ and 
because 717- is an epimorphism, v + v = v s + v £ , where v s G S and v £ <E £. 
Hence, v — v s = v £ — v , where v — v s G S and v £ — v G £. Since £±5, we 
arrive at a contradiction. 

(v) Similarly, we show that 7rJ^(w) = (w + J 7 ) fl J\f, where T = (S fl 
T) x n S = S', and dim J 7 = dim£. 

(vi) Take a cube of radius r in S which is the direct product of a cube 
of radius r in T C\S and a cube of radius r in iS'. Then its image equals the 
direct product of a cube of radius r in T fl S and a subset of T. Hence, it 
suffices to calculate the measure distortion of -k^- 1 on S' . Similarly, we infer 
that it suffices to calculate the measure distortion of irj^ 1 on £'. 

(vii) Consider the orthogonal mapping 6 that is a rotation about span{A/"fl 
£, T fl S} such that 

• G(F) = £ 

• e((A r n £) x n £) = (T n S) 1 - n 5. 

Note that if we choose an orthonormal basis in R N in the following order: 
span{A/" fl£,Tn5}, spanjj 7 }, span{(A/" fl £) _L fl £} then the matrix of 6 
looks like 

/ E 
0=0 E 
V £ 



where 



/ 










1 



o \ 





10 



\ 1 / 

Consequently, 7r r 1 \ s/ (B(0,r)nS') = o vr^ 1 \ C >(B(0, r)n£') for anyr > 0, 
and the measure distortions of these mappings coincide. 

(viii) Thus, the measure distortions of and ttj- coincide. 

Since both the Riemannian differential and foe-differential of 6 X at x equal 
identity, by the chain rule the corresponding "determinants" coincide: 

V(D<p(x)) = P(£ty(0)) and V(D<p(x)) = V(Dijj(0)). 
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Thus, the measure distortion under 7r r equals 

V(Dip(x)) _ V(Dif(x)) 
V(Dijj(x)) ~ V(Dip(x))' 

Step III. Verify that ^~ 1 (t)nBox 2 (0, r) is a subset of an o(r)-neighborhood 
of ker Dip (0) fl Box2(0,r(l + o(l))) with respect to d 2E . To this end, we 
prove that <9(?/>~ 1 (t) fl Box 2 (0,r)) is a subset of the o(r)-neighborhood of 
<9(ker Dip(0) flBox 2 (0,r(l + o(l)))) with respect to d 2 , where o(l) is uniform 
in x = 9 X {0), x G U m M. 

Indeed, consider y G fl Box 2 (0, r)). We can represent this point 

as y = y Slp + y kcv 3r where Vp* e (kerD^(0)) ± and y kci3i) G kerZ^(O). 
Then the definition of the /ic-differentiability implies that 

o(r) = di {o \5m(y),m) = 4 (0) w(o)(^)^(2/)) 

= 4 (0) W(o)(^),V(o)). 

We have 4 (0) (^(0), 5^(0)(^)) = o(r), and ^(0,y 5l/) ) = o(r), where o(-) 
is uniform in y, x = 9 X (0), x G U <<= M. The inequality 

d° E (0, « + v) < d° 2E (0, u) + cd° 2E (0, v) (3.9) 

implies that 

4E(^y kcr D^< d 2E(^y) + cd° 2E (0,y^)=r + o(r), 

and hence, d° 2E (0, y kerSl/) ) = r(l + o(l)) and d°(0, y kerSv ,) = r(l + o(l)). Thus, 

d(^ l {t) fl Box 2 (0,r)) is a subset of the o(r)-neighborhood of d(ker Dijj(0) fl 
Box 2 (0,r(l + o(l)))) with respect to d 2 - 

Step IV. Similarly considering the linear mapping L(y) = Dip(0)y, we 
infer that <9(ker Dip(0) D Box 2 (0,r)) is a subset of the o(r)-neighborhood of 
<9(ker Dtp(0) fl Box 2 (0,r(l + o(l)))) with respect to d 2B . Indeed, we have 
ker£ty(0) = £ _1 (0) and DL(0) = Dip(0), and L is /ic-differentiable at 
since Z>0(O) is horizontal, and -DL(O) = Dip(0). 

Step V. Since there exists a bijective linear mapping from kerZ>0(O) fl 
Box 2 (0,r) to kerDip(0) n Box 2 (0,r(l + o(l))), by Step IV we infer that 
<9(ker _D-^(0)nBox 2 (0, r)) is a subset of the o(r)-neighborhood of <9(ker Z>0(O)n 
Box 2 (0,r(l + o(l)))). 

Hence, c^ -1 ^) fl Box 2 (0,r)) is a subset of the o(r)-neighborhood of 
<9(ker Dip(0) n Box 2 (0,r(l + o(l)))) with respect to di^. 
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Step VI. In Steps VII - IX we explain that 

d(kexDip(0) nBox 2 (0,r(l + o(l)))) 

is a subset of the o(r)-neighborhood of c^ -1 ^) fl Box 2 (0,r)) with respect 
to d 2E . Moreover, similar arguments imply that the same is true regarding 
the sets ker Dip (0) n Box 2 (0,r(l + o(l))) and ip~ x {t) H Box 2 (0,r). Indeed, 
it suffices to recall that o(-) in the previous steps are uniform in r > and 

x = e x (o). 

Step VII. In order to justify the result of Step VI, we construct a one-to- 
one mapping from <9(?/> _1 (£)nBox 2 (a;, r)) to a subset of the o(r)-neighborhood 
of <9(ker Dip(0) flBox 2 (x, r(l + o(l)))) (with respect to <i 2£; ) lying in the plane 
kerD^(O). 

First of all, observe that at the regular point we have ker Z)^n(ker Dip) 1 - = 
{0} by Step II, substep (iii). 

From this we deduce the following assumption: without loss of generality 
we may assume that Dip(z)(v) > (3, where z G U D Box 2 (0,r), (3 > 0, for 
v G (ker Dip (0)) 1 - with \v\ = I. 

This implies that in a neighborhood of a regular point we can expand 
each point y uniquely as y = ykerD^ + VQipi where ykerDV G ker Dip and 
ygw, e (kerD^)- 1 . Suppose that y G ip~ l {i) fl <9Box 2 (0,r) and consider the 
mapping £ : y i-)- y k er £></>• 

Verify that <i 2 (0, £(y)) = r(l + o(l)). To this end, we prove that <i 2 (0, y — 
i(y)) = o(r) and then apply (13. 9 p for d 2E . 

Represent y as y = y kC rD^ + 2^ = v ker dip + v d^ + where V^d^, = 
v ker5^ + v 5r Since d 2E( G ,y) = r > Ste P 111 implies that dl E (0,v keiSf ) = 
r + o(r) and c? 2£; (0, vq , + yg^) = o(r). Moreover, by Step IV we have 

<^b(0, = (1 + o(l))d° 2E (0, v keiBip ), d° 2E (0, v Dil ) = o(l)d° 2E (0, v kerBi> ), 

and taking into account the generalized triangle inequality for d 2E , we con- 
clude that d 2E (0, y Bl p) = o(r). Finally, we have d 2 (0,y — £(y)) = o(r). 

Step VIII. On this step, we show that £ is a bijective mapping. It is easy 
to see that £ is injective. Indeed, it follows from Dip((ker Dip) ± fl S 7 ^ 1 ) > 
/3>0. 

To verify that it is also surjective, we show that the mapping 77 : ip^ 1 fl 
Box 2 (0, 2r) — > ker Dip(0) is bi-Lipschitz. Here, 77 assigns to each y G ip^it) fl 
Box 2 (0,r), the vector y kerD ip, where y = y ker D^ + yg^, ykcrD^p G kerDV, and 

y^ G (kerfh/0" 1 . It is clear that »7|^-i(t)naBax3(o,r) = £■ 

Since 77 is a projection, it is Lipschitz with the Lipschitz constant being 
equal to 1. Verify that |^(yi) — 77(2/2) | > ^|yi — 2/2 1 for some K > for all 
2/1,2/2 <E V _1 (*)nBox 2 (0,2r). 



36 



Suppose that it is not so; thus, for every e > there exist y 1: y 2 € ip 1 (t)r\ 
Box 2 (0,2r) such that ^(j/i) — ^(2/2) I < e\yi - y 2 \. Since ^ = y ikcrDi , + yi^, 
Hiker Di, e kerD^ and 2/jg^ G (ker Dip)- 1 -, % = 1,2, we infer that 

\v(yi) - v(y2)\ < s\yi - 2/2I < eWs/i) - 17(2/2) I + e|?/is^ - 2/25^1, 

and 

- 2/25^1 ^ — ^ — 1^(2/1) - 17(2/2) I- 

Since kerZ>0(O) fl (ker /^(O))" 1 = 0, the norm I2/1 — 2/2^ = \v(yi) ~ ^7 (2/2) | + 
\yiDiit ~ y^D^l * s equivalent to the Euclidean norm. Furthermore, 

2H2/1) = ^(2/2) + Dip(y 2 )(y 1 - 2/2) + o(|yi - y 2 \), 

and Dip(y 2 )(y 1 — y 2 ) = o(\yi — y 2 \). Since 77 is a linear mapping, we have 

2/15^ -yzBi, = (yi - y^D^ and ^(j/O - vfa) = v(yi - 2/2). 

On the one hand, by the definition of i] we have 

^(2/2X2/1 - 2/2) = £ty(0)(j/i - 2/2) + oflyi - 2/2I) 

= £ty(0)((yi- 2/2)^) + o(|2/i-2/ 2 |), 

and D^(0)((j/i - 2/2)5^) = o(|j/i - 2/2 1 ) - On the other hand, 



-2/25^1 > ^— +1 J JI2/1-2/2U > L(i-{—— +ij JI2/1-2/2I, 

where L > depends only on U. 

These relations lead to a contradiction with 

DV(0)((ker5^(0)) ± ) > /3 > 0. 

Thus, i] is a bi-Lipschitz mapping, and each 2/ £ ^^(t) n<9Box 2 (0,r) has 
a unique preimage. Thus, £ is bijective. 

Step IX. Step VIII implies that £ is invertible, and the previous steps 
imply that the ^-distortion of is (l+o(l)) with respect to 0. The estimate 
for c? 2 (0, 2/ — ?7 1 (2/) ) is proved in Step VII. Observe that o(l) depends on the 
convergence of o(l) to in the equality 

4 (0) (^(0)2/,^(2/)) = o(l)-rf 2 (0,2/). 

Thus, o(l) is uniform in y. Consequently, V' _1 (t) n<9Box 2 (0,r) is a subset of 
the o(r)-neighborhood of ker .D-?/>(0)n<9Box 2 (0, r(l+o(l))) with respect to d 2 , 
and conversely, ker Dip(0) n<9Box 2 (0,r) is a subset of the o(r)-neighborhood 
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of ip^it) n <9Box 2 (0, r) with respect to d® E (see Step VII). The same is true 
regarding the sets ip^it) fl Box 2 (0,r) and ker-D^O) D Box 2 (0,r(l + o(l))). 

Note that the extension rj of £ onto ^ _1 (/!;)nBox 2 (0, r) has the ^-distortion 
with respect to equal to 1 + o(l) as well. In Step VIII, we prove that this 
extension is bijective. 

Step X. On this step, we show that the classical, that is, "Riemannian" 
distortion of r\ (and of £) is also 1 + o(l). 

Indeed, since kerD^(O) n (ker D^(O))- 1 = {0} and L>^((ker 5^(0)) ± ) > 
/3 > on [/, it is easy to see that p(0,yg^,) = o(p(0, y)). Hence, p(0,yg^) = 
o(p(0, yk er D?/))), and the p-distortion of both 77 and r/ _1 equals 1 + o(l). 

Since o(-) are uniform in the definitions of Riemannian differentiabil- 
ity and /ic-differentiability, o(-) is uniform (both in the case of d?r and p- 
distortion). 

Thus, claim I is proved. 

Step XI. Using the results of all previous steps, we can consider the 
mapping 7]~ l (that is, the extension of £ _1 ) from ker Dip n Box 2 (0,r(l + 
o(l))) to ip~ l {t) flBox 2 (0,r) (see Steps VII IX). The established properties 
of r) imply that the measure of ij) (t) fl Box 2 (0,r) equals that of ker Dip fl 
Box 2 (0, r) up to a factor of 1 + o(l), that is, it equals (7(1 + o{l))r u ~ v , where 
C is obtained above. 

Indeed, rj is a bi-Lipschitz C^-mapping since it is a projection of a C 1 - 
surface in a nontangent direction. Consequently, the inverse mapping r]~ l 
is also a bi-Lipschitz C 1 -mapping. Moreover, since it is differentiable, it is 
also metrically differentiable, and its metric differential equals 1 + o(l) on 
every direction. The C 1 -mapping r]~ l satisfies 1-L N ~ N (r]~ l (A)) = J7"(r/ _1 ,0) ■ 
n N ~ N (A) ■ (1 + o(l)), where A = 77(V> -1 (t) n Box 2 (0,r)). We can calculate 
the Jacobian of rj^ 1 via its metric differential [Kil IKml} IKm4| IKm3| IKm5j : 

J( V -\0) = <r N -s[ J (l + o(l)f- N dH N -^)]'' 1 = (l + o(l)). 

§JV-iV~l 

Thus, H N - N (ji^^A)) = H (A) (l+o(l)). By taking into account a remark 
in Step V substep (ii) of the proof of Theorem 13. 71 concerning regular points, 

M 

we infer H N ~ N (A) = U^n^ • gigf^ 1 + 

Step XII. Consider the restriction of 9 X to ip~ l (t) fl Box 2 (0,r). Observe 
that eJlijj- 1 ^) nBox 2 (0,r)) = y~ l {t) n Box 2 (x, r). 

The "H^^-measure distortion under this mapping equals X > (<7|ker.D<p(a;))> 

and CRiem = T>(g\keiD<p(x))- 

Claim II is proved. 
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All o(-) in both claims I and II are uniform in the radius r and x = 6 X (0), 

x g U m M. 

The proof of the theorem is complete. □ 

Definition 3.12. The (spherical) Hausdorff measure of a set A C <fi~ 1 (t) 
(constructed with respect to a sub-Riemannian (quasi)metric d and sub-Rie- 
mannian balls Bd in d) equals 

W-\A) = ^liminfj^rr^ : \J B d (x h r<) D A,x t e A,r t < 6,1 e n}. 

Property 3.13 ( |KmVod] ). The quasimetric d 2 and the metric <i cc (see Def- 
initions [223 and [2UI]) are locally equivalent. 

Lemma 3.14. Given a set A C y? -1 ^) of T-L N ~ N -measure zero consisting 
of regular points and e > 0, there exists a covering of A by the "balls" 
{Box 2 (xj, Tj) PI <£> -1 (t)}j e N; where Xj G <^ _:L (i) ; j G N, the sum of whose 
H N ~ N -measures is less than e. 

Proof. Fix e > 0. Represent A as the union of some subsets Ai C A, I G N, 
lying at positive distances from the set \ H <£> _1 (t). Without loss of general- 
ity we consider a set A; instead of A. Suppose that Ai satisfies the stated 
conditions. Observe that for this set there exists a collection of "balls" 
{B(xj,rj) fl where Xj G <£> -1 (t) and B(xj,Tj) are Riemannian 

balls, j G N, the sum of whose "H^^-measures is less than e, where e is 
determined by e (we specify the exact expression for e below). Fix j G N. 
For B(xj,Tj) fl (p~ l {t), there exists a collection 

{B cc (x, r) n ^(t) :ieA,n ¥> -1 (t), 5 cc (x, r) n yT^t) C B{x j ,r j ), r > 0}. 

By the 5r-covering lemma (see [Fe2j for instance) there exists a countable 
family of disjoint "balls" {B cc (xj k , Tj k ) fl ip~ l {t)} such that 

fceN 

(Here we use the metric <i cc since the 5r-covering lemma is established for 
metrics, and we cannot say for sure whether it holds for quasimetrics) . Since 
d 2 and d cc are locally equivalent, there exist some constants < Ci, C 2 < oo 
such that Box 2 (x, C\r) C B cc (x, r) C Box 2 (x, C 2 r) for sufficiently small r > 
and x in some sufficiently small neighborhood. Consequently, there exists a 
disjoint collection {Box 2 (x-, fc , Cirj k )} such that 

[J Box 2 (xj k , 5C 2 rj h ) D Ai H B( Xj, Tj). 

fceN 
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This implies that 

]>>^(Box 2 (^ fc ,5C< 2 r,J) < C(5, C lf C 2 ) £«"-*(Bcx 2 (s ilb , C x r,J) 

fcGN fceN 

< (7(5, Ci, C 2 )n N - N (B( Xj , r,-) n 

Then {Box 2 (x., fe , 5C 2 r, fc )}j fc6N is a reqired collection, and the sum of the 
-^7v-iv_ measures f these "balls" is at most C(5,Ci,C 2 )e < e. □ 

Corollary 3.15. For each % N ~ N -measure zero subset o/<^ _1 (t) consisting 
of regular points we have T-i v ~ u = 0. 

The proof follows directly from the definition of an 'H'^-negligible set. 

Lemma 3.16. Consider a regular point x G M. Then given a sufficiently 
small r > with Box 2 (a;,r) fl <£> _1 (<£>(V)) fl % — an d £ > 0, there exists a 
covering o/Box 2 (x,r) fl by the sets Box 2 (xj,rj) fl (p~ l (<p(x)) with 

Xi G and Box 2 (xj,rj) C Box 2 (x,r) ; i G N, such that 

Y,n N -"(BoMxi,ri) n^r 1 ^))) < H Ar -^(Box 2 (x, r) n ^\^{x))) + e. 

Proof. Without loss of generality we assume that r > satisfies 
dist(Box.2(x,r)r\(p~ 1 ((p(x)),x) > 0. 

Since the measure % N ~ N is doubling on Box 2 (?/, s) fl (/? _1 (<£>(:r)) for y G 
(p~ l ((p(x)) and sufficiently small s, in view of the Vitali covering theorem 
there exists a collection of "balls" {Box 2 (x ij , r^.) fl <£> -1 ((/?(:r))}j eN with G 
<p -1 (<p(x)) and Box 2 (x ij , r^) n <^ _1 (^(a;)) C Box 2 (x,r) n <p _1 (<p(:r)), j G N, 
such that 

^-^(Box 2 (x,r) n^ 1 ^))) = ^^ JV - J? (Box 2 (x i .,r ii ) n<p-\<p{x))). 

By Lemma [3.141 the ^^"^-negligible set 

Box 2 (x,r) n ip~ x {<p(x)) \ (J Box 2 (a;i j ,r ij .) n ^(^(a;)), 

jen 

admits a covering by the collection {Box 2 (xj fc , r ik ) fl <£> -1 (<£>(V))} of "balls" 
the sum of whose "H^'^-measures is less than the given e > 0. The proof is 
complete. □ 
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Recall that D til /j l2 (y) stands for the derivative of a measure /i 2 with respect 
to a measure \i\ at y: 

D^n 2 {y) = hm— — -. 

Theorem 3.17 (The Lebesgue Differentiation of Measures on Level Sets). 

The Hausdorff measure % v ~ v of Box 2 (x,r) fl (p^^ip^x)), where x is a reg- 
ular point, and dist(Box 2 (x, r) fl (^j -1 (</?(#)), x) > 0, asymptotically equals 

n v - f> (Box 2 (x, r) n f-\f(x))) = Uv-pr^il + o(l)). 
The derivative -Da,]v-iv U v ~ v {x) equals 

1 Uu-v _ X>(Dy?(g)) 

^ (Piker D»)(i 

11 ^rik-rik 
k=l 

Proof. Consider the intersection of the sub-Riemannian ball of a sufficiently 
small radius r > centered at x and the level set ip~ 1 ((p(x)). Fix 5 > and 
some covering {Box 2 (?/j, rj)}j e N of this intersection by sub-Riemannian balls 
as in the definition of the set function H$~ u ■ Recall that for a (quasi) metric 
d we have 

n u s -°(A) = w„_p inf {£)rr e = [jB^n) DA,x t e A,n < 6,i G n} 
Then, setting 



a(w) = lim — 

H N - N (Box 2 (y, r) n ^^(y))) 

by Theorem 13.111 we have 

= + %^)]^^(Box 2 (y i ,r i ) nyr 1 ^*))) 

+ A(x,r)]^H 7V -^(Box 2 (^,r,)n^- 1 (^(a;))), 



L a l x J 



where <5(?/j, r*) — > as r 4 — > uniformly in z e N, and A(x, r) — >■ as r — >■ 0, 
since .D<£> and .Dy? are continuous, and by Theorem 13.111 as well. For the 
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fixed x and e > 0, we can choose a sufficiently small radius r > to satisfy 
\A(x,r)\ < e ■ a(x). Consequently, 



(1 - e)a{x) J2u N -*(Bo X2 ( yi , n) n <p~\<p(x))) 

igN 

< u^J2 r ^ ^ t 1 + ^^^(60x2^,^) n ^\^{x))). 

Note that for each 5 > if the covering in the definition of 7ig~" is nearly 
"optimal" (i. e., r i~ u * s near ly minimal), then % JV ~ JV (Box 2 (?/j, r,) fl 

is also nearly minimal. Since the T^'^-measure is countably 
additive, by Lemma 13.161 the infimum of the values of the last sum equals 
n N ~"(Box 2 (x,r)). 
Hence, 

n u -%Box 2 (x, r) n <p- l (<p(x))) = uw^(l + o(l)), 
where o(l) — > as r — >■ 0, and 

D vN - i iU v -*(x) = a(x) = lim = (3.10) 

n ^o-H^- Ar (Box 2 (x,r)n^- 1 (^(x))) 



1 u v - s V{Dip(x)) 



^ (Piker D<p(x)) 

11 ^rik-rik 



(3.11) 



fe=i 

The proof is complete. □ 

The last result motivates the next definition. 
Definition 3.18. The sub-Riemannian coarea factor equals 

J~ (<p, x )=V{D<p(x))-— — — . 

v N TT 

11 "'rifc— n* 
k=l 

Remark 3.19. By (ETTOj) . we have 

jfV, x) = J^, x) . ■ ■ ^ 

X>(g>(^(x))) ^ jv 

since _ _ 

V{g{x)) 
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Corollary 3.20. Theorem\3 '. 1% Corollary \3.15[ and the Lebesgue differenti- 



ation theorem (see, for instance, |VodU] ) imply that every measurable subset 
of a regular set A C (D D ip^it)), t 6 M, satisfies 

H V -\A) = J D nN _ f} n u - D (x)dn N "^(x). 

A 

Theorem 3.21 (Local coarea formula for the set D). If (p G C^IM^M) then 
Jl R ( V ,x)dH v (x) = J dU\t) J dU v -\u). 

Proof. The Riemannian coarea formula [Fe2j yields 

J£ R (<p,x) 



j J§ R (<p,x)dU l '(x) = j * * J dU N {x) 



du N (t) I J l R ^z u) du N -»{u). 



J R (<p,u)V(g(u)) 

M <p _ (*)no 



Calculating the measure derivatives in M, we obtain 

J J^,u)V(g(u)) 

p- 1 (t)nB 



M ^ _1 (*)nB 

Finally, taking Remark 13 . 1 9 1 and Corollary 13.201 into account, we infer that 



J R (<p,u)V(g(u)) 

M (p-Ht)nB 

= y d^^t) y du v -\u). 

The proof is complete. □ 
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4 The Characteristic Set 



The goal of this section is to prove 
Theorem 4.1 (see also |Km7j ). If Assumption is fulfilled then 

u v -\ X t) = o 

for almost alltEM. (with respect to both T-L N and ~H U ). 

Recall that we denote ip H x D Y Xt f° r t G M. 

We prove the results of this section under the following assumptions. 

Assumption 4.2. Assume that the basis vector fields {A^}^ in the preim- 
age are of class C M+1 , and {Xj}^ =1 in the image are of class C 1,a , a > 0. 

Moreover, assume that <p G C A/+2 (M, M), where M is the depth of M (see 
Definition ED}. 

Lemma 4.3. Let M be a Carnot manifold with C 1,a -smooth basis vector 
fields, a > 0. From the basis vector fields {X{\f =l construct linearly inde- 

dim H dEg x . 

pendent vector fields {Yj}^ =1 as Yj(x) = Yl a ji(x)Xi, where dji(x) are 

8=1 

C 1,UJ -smooth functions for all i and j , u> > 0. Then the quasimetric of type 
g?2 constructed with respect to the new basis {Yj}j =1 is locally equivalent to 
the initial quasimetric d 2 . 

Proof. Denote the quasimetric of type d 2 and the Carnot-Caratheodory met- 
ric constructed with respect to the basis {Yj}^ =1 by d\ and o% c respectively. 
The Ball-Box Theorem of |KmVodl IKm9j implies that d 2 and d cc are lo- 
cally equivalent. Moreover, since the lengths of a curve calculated with re- 
spect to different bases are bi-Lipschitz equivalent, the metrics d\ and d Y cc 
are bi-Lipschitz equivalent as well. Considering that d cc and d\ c are locally 
equivalent, we infer that so are d 2 and d\. □ 

Proof of Theorem \4-l\ Observe that the condition % N (x) — W(x) = an( i 
\Fe2\ 2.10.25] immediately imply the result. 
Assume now that 7i N (x) > 0- 

Step I. Divide x m to finitely many sets, each of which is defined by the 
structure of ker Dep. In other words, represent x as 

X = Xpx—PMI 
pi,—,PM 
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where 



Xpi-pm = \ x e X '■ dim(Hi(x) n ker Dip(x)) = p i: i = 1, ...,Mj. 

The assumption T-t N (x) > implies the existence of a collection {pi, . . . ,pm} 
with H N ( Xpi ... PM ) > since ft"( x ) < £ 

Step II. Fix a collection . . . ,Pm} with this property, the correspond- 
ing set Xpi-pm: an d x ^ Xpi-pm (which is not necessarily a density point). 
The purpose of Steps II-V is to prove that 

• the set Xpi-pm * s measurable; 

• in the coordinates of the first kind the order of tangency between a 
level set and its tangent plane is o(r M ) at almost every x G Xpi-pm 
(with respect to both % N and W), 

these facts that are independent of the basis transformations described in 
Lemma 14.31 on sufficiently small U (s M. Thus, without loss of generality we 
may assume that at this chosen point x the differential Dip(x) is degenerate 
on basis vector fields X^x), . . . , ^(^)- Consider a compact neighbor- 
hood f/iMat whose points y the differential Dip(y) does not vanish on the 
same vector fields as at x. 

Step III. On Step IV, we are going to prove that % N (Xpi~pm) ^ s mea_ 
surable. To this end, we now construct a new auxiliary basis {Yj}jLi in the 
same way as in Lemma [4.31 This construction consists of several steps. 

(i) The purpose of this substep is to obtain a new horizontal basis the 
nonzero images of whose vectors under Dip at X P i...p M are independent. 

The definition of Xpi-pm anc ^ ^ ne structures of Dip and Dip (see ( 12 .5p and 
(J23D) imply that 

dim(ker Dip(x) fl Hi(x)) = dim(ker Dip(x) D Hi(x)) = pi, 

where p\ > rii — hi (see Proposition I2.21[ claim I (a)). Then, the images of 
the remaining n\ — p\ horizontal vectors are linearly independent on ip(U). 

Next, we consider the matrix of the classical differential Dip written with 
respect to the basis vector fields Xi, . . . , Xn in the preimage and X\, . . . , X^ 
in the image. 

Consider the hi x rii block of the matrix of Dip that corresponds to 
the horizontal subspaces. Identify its ni columns as vectors in K m . Then, 
at x we choose from this collection rii — pi linearly independent vectors 
Wx(x), . . . , w ni - Pl (x). At a point y G U we denote the corresponding column 
vector as Wi(y), . . . , w ni - Pl (y). The remaining p\ vectors V\(y), . . . , v Pl (y) 
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at y G Xpi...p Ml ^ U are contained in the linear span of the first ones. 
Consequently, for the set X Vx ... VM there exists a transformation 0™1- P1 (y) 
of M ni depending C M+1 -smoothly on y G U such that the images of all 
Wi(y), . . . , w ni - pi (y) belong to x o ni- ( ni-Pl ). To see that this map- 

ping is C M+1 -smooth, it suffices to write the coordinates of the linearly 
independent vectors {wi(y), . . . , w ni - Pl (y), e h , . . . , e ini _ (ni _ pi) } in the stan- 
dard basis as a matrix and then to calculate the inverse matrix. (Here 
{e^, . . . , ei ( } } are the vectors of the standard basis in IR ni linearly 
independent together with wi(x), . . . , w ni - pi (x); without loss of generality 
we may assume that the neighborhood U is sufficiently small for the vectors 
{e h , . . . ,e ifti _ (ni _ pi) } to be independent together with wi(x), . . . , w ni - Pl (x) 
for all y G U.) The resulting matrix is the matrix of this transformation. In 
other words, the matrix 

(OZ_ pl (y)[My)l Ol\_ pi {y)[w ni - Pl {y)\) 

has exactly n\ — {n x — pi) vanishing columns. Since we can apply to these 
vectors an orthogonal transformation independent of y G U without loss of 
generality we may assume that these vanishing columns have indices from 
ni — pi + 1 to n\ . 

By the linear dependence, at the points of Xv\-vm ^ ne images of the col- 
umn vectors vi(y), . . . ,v Pl (y) also belong to M ni ~ Pl . Represent these images 
(at an arbitrary point y G U) as 

ol\- P M[v k {v)\ = H{y)^l{y)) T , 

where the dimension of oo\(y) equals ni—pi, and the dimension of (jjf.(y) equals 
hi — {rii — pi). Since 0™J_ pi (y) G C M (U), where U is a compact neighbor- 
hood of the origin, X; G C M+1 , % — 1, . . . , JV, it follows that io\{y) depends 
C M -smoothly on n\ — p\ columns. Hence, there exist smooth coefficients 

ni-pi 

6 M (y),...,6 fcj7ll _ Pl (y) with u\(y) = b k ,i{y)wi{y) . Then if the columns 

1=1 

w^y),..., w ni - pi (y) correspond to the vectors X h (y), . . . , X ini _ pi (y), and a 
column Vk(y) corresponds to Xj k (y), then upon replacing Xj k (y) by 

m-pi 

x jk (y)- E Ki(y)x k (y) (4.1) 
i=i 

we find that the part of the new column corresponding to u)\{y) equals 
zero, where k — 1, . . . ,p 1 . Next, since at the points of Xpi...p M ^ ne columns 
Mv), ■ ■ -i v Pi{y) belong to M" 1_Pl we have u 2 k {y) = on x Pl ... PMl - A PP l Y the 
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transformation (14. Xp for all k — 1, . . . ,p±. In this case, each new column 
v k (y) vanishes at y G Xpi- PM i k = l,...,p 1 . 

Thus, we obtain a new horizontal basis the nonzero images of whose 
vectors under Dip at Xpi-p M are independent. 

(ii) Next, we apply similar arguments to the subspace H 2 using the fact 
that we have already constructed a basis for Hi. Namely, at y G Xpx—pui 
in each of p 2 — Pi vector fields, we delete the parts depending on n 2 — p 2 
linearly independent vector fields. This yields a new basis for H2 the nonzero 
images of whose nonhorizontal vector fields independent and nonhorizontal 
(belonging to H 2 ). 

Continuing similarly, we obtain new bases for all H^, k = 3, . . . , M. De- 
note these vector fields by {Yj}^ =1 . 

Step IV. By step III, the intersection of Xpi...p M w hh the closure of every 
compact neighborhood W <s U is closed. Indeed, take {yi} —> y as I —> 00, 
where yi G Xpi-pm ^ W. Since Dip(yi) at each point yi vanishes on a certain 
collection of basis vector fields by the definition of Xpi-pm 1 ^ ne continuity of 
the differential implies that D(p(y) vanishes on the same vector fields and is 
non-degenerate on the rest since W <<= U. Thus, the set Xpi-pm * s measurable. 

The measurability of Xpi-pm implies that almost all its points are density 
points (both with respect to T-L N with balls in the Riemannian metric and 
W with Box 2 -balls). 

Step V. Estimate now the "H^'^-measure of the intersection of the level 
set passing through a density point x G Xpi-pm °f Xpi-pm 1 anc ^ ^he ball 
Box 2 (x,r). To this end, we use in the image and preimage the normal co- 
ordinates with respect to the points (p(x) and x, and estimate the order 
of tangency between this level set and the tangent plane. We denote the 
resulting composition o ip o 6 X by ip. 

N-N N-N 

Evaluate ip at an arbitrary point y = YI UZjXfy = S ^i e jn where 

i=l i=l 

Zj t = D8~ l [Yj.] are the basis vector fields tangent to the level set ^> -1 ( , 0(O)) 
at zero. We have |NSW| 

M N-N 



^) = E[t[ (E ^)V](0) + o{\y\ M ) 

3=1 J ' i=l 

M 1 N-N 

= E[7f(E^)V](°) + °(H M )' 
3=2 j ~ i=i 

where o(-) is uniform on U. Since is a density point of zeroes of DipZ^, the 
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result of the action of every differential operator on it vanishes at as well: 

N—N 



(X> z *) iKo) = o. 



Similar statements hold for almost all (with respect to both TL and %") 
points of Xpi—pm' Consequently, is a density point of the zeroes of the 
function 

N-N 2 

Then we obtain 

M , N-N 



^) = E[7t(E^)^](°) + °(H m ) 

j=3 J ' i=l 

M N-N 

-Eb(E^)^] (0)+o(l2/|M) 

j=4 J ' i=l 

,N-N \ 3 

since £ UZ k V(0) = 0. 
v i=i ' 

Applying similar arguments, we arrive at ^(y) = o(\y\ AI ). This implies 
that the tangent plane approximates the level set up to o{\y\ M ). Indeed, 
take y in the same level set as 0, and the orthogonal projection y' of y along 
(ker DipifS)) 1 - to the tangent plane. Then 

<%) = Hy') + o(\y\ M ) = Hy') + Di;(y')(y-y')+o(\y-y'\) = ^y') + C\y-y'\, 

and since |C| > C > uniformly on U, we have \y — y'\ = o{\y\ M ). Observe 
that here o(-) is uniform on U. Consequently, Theorem 13.71 yields 

n N-N^-i^ Q ^ n BoX2 ( , r )) = Lr v ~ v \ 
where 0<L 1 <L<L 2 <oo uniformly on U, and 

H N -"{ip-\ip{x)) n Box 2 (x,r)) = ar"-^, (4.2) 

z^o > ^ Here the number v is defined by the set Xpi...p M 'i n d e P en dently of x. 
Since the order of tangency is independent of the basis, the box in (14.21) is 
taken in the initial quasimetric di which is constructed with respect to the 
initial basis {X^\f =l . Moreover, the relation 

n N -"(<p- 1 ( V (y))nBo X2 (x,r)) = C y r^ 
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with < K\ < C y < K 2 < oo holds for all points y G X Pl ... PM \ £ Xpi p for 
some set S y with H N '(S y ) = ^(E y ) = 0. 
Step VI. Suppose that zeM. Estimate the measure 

(z) n (x P1 ... PM \ s Xpi ... pm ) n wo 

for W U. Since d 2 and rf cc are locally equivalent, the conditions of 5r- 
covering lemma are fulfilled for the balls centered at the points of Xpi...p M \ 
£ Xp p . Consider a covering [j Box 2 (xj,rj) of the set 

^ 1 Wn(x P1 ... PM \s Xpi ... PM )n^ 

such that the balls {Box 2 (xj, Tj/5/)} ie N are disjoint. Here the number / de- 
pends on the constants in the equivalence of <i 2 and d cc . Then we have 

U v r^p-\z) n (x Pl ... PM \ s Xpi ... pm ) n w) 

< c^'-^r^-^ < c^-^ft"-^- 1 ^)) n 60x2(0;,, n)) 

< C 2 ^- v Y^n N ~ S &~\v( x )) n 60x2(^,^/5/)) < 8 V °~ V L -> 

as 5 — )■ since L < H N ~ N ((p~ 1 (z) fl (x Pl ... PM \ £ Xpi p ) fl 14 7 ) is independent 
of 5 > 0. Thus, the intersection of every level set with Xpi-pu \ ^Xpi...j> nas 
% v ~ v -measure zero. Consequently, 

^-%- 1 (^)n( x \s)) = o 

where E = |J £ Xpi ... pm and ^(S) = ft"(£) = 0. Now the theorem 

Pl—PM 

follows from [Ml 2.10.25]. □ 

5 The Degeneration Set 

Theorem 5.1. For TL" '-almost all t G M we have 

u v - p (<p-\t)nz) = 0. 

Proo/. Step I. We have Z = Z 1 U Z 2 , where VP{^{Z X )) = U s {ip{Z$) = 0, 
and 'H N - N {(f- 1 {t) nZ 2 ) = for all £ G M (see, for instance, [Km6]). 
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Step II. Note that % N (<^(Z 2 )) < oo since otherwise we can represent it 

iv-i 

as a countable union of sets of finite measure. Next, decompose Z 2 as [J Ci, 

i=0 

where rank Dip{x) = i for x e Cj. 

Fix < i < N — 1. Without loss of generality we may assume that Cj is 
a compact set. Consider a compact set Cj C C such that <p(Ci) is a compact 
subset of the set of density points of <p(Ci) and l-L N ((p(Ci) \ <f(Ci)) < e for 
a fixed arbitrary e > 0. Consider the mapping 

= xjv_ i+ i, . . . , xat) : — >■ IR^ x R N ~\ 

Observe that £j(C) is an iV-rectifiable set. Our assumption implies that £j 
is locally bi-Lipschitz on Cj. Consequently, % N ~ N '(£i(<^ _1 (£) PI Cj)) = for 
all tel. Consider the mapping 7r : £(Cj) — > R defined as 

Apply the classical coarea formula to 7r, Cj and C: 
J J(ir,x)dn N (x)> J XtfaWJfax)^"^) 

= / «*H*(f) J Xw^dU^iu). (5.1) 

^(Cj) 7r-i(i)nc(C0 

By the definition of Cj, we find that J7"(7r, x) = 1 for x G £(Cj). Furthermore, 

^-^(tt- 1 ^) n £(C)) = W^Ce*^- 1 ^) n C)) = o. 

It follows from (JSH]) that H N (£(&i)) = 0. Finally, have 

?^(Cj) = ?T(C) = 

because £j is bi-Lipschitz. Since e > is arbitrary, there exists a set Ej of 
measure zero in ip(Ci) such that 

u N (v~\y(Ci) \ Ej) n c^ = H v {<p-\<p{Ci) \ Ej) n d) = o. 

Therefore, W[ U y?"%(Cj) \ Ej) fl Cj) = 0. 



50 



Step III. Since cp is contact, it is locally Lipschitz with respect to every 
sub-Riemannian metric. Reasoning as in |Fe2[ Theorem 2.10.25], we obtain 



0<y dH*(t) J dU v -\u) 

M JV-i 

rH*)n U v- 1 (v(c , l )\s i )nCi 

i=0 



JV-1 



< C(Lip(p), u, v)W ( |J y-\y{c % ) \ so n Ci) = 0. 

i=0 



JV-1 



Furthermore, the T^-measure zero set |J Ej is included into an T^-measure 

i=0 

zero Borel set E. Then, its preimage <^ _1 (E) is also a Borel set, and the set 
<^ _1 (E) fl Z 2 is measurable. By the definition, 

H"(<p(<p-\E) n z 2 )) = ^(^-'(E) n z 2 )) = o. 

combining the results, we deduce that 



dU v {t) J dU v - v {u) = 0, 
and thus we complete the proof. □ 

6 The Coarea Formula 

Theorem 6.1 (Coarea Formula). If a contact mapping (p : M — > M satisfies 
Assumption \2.2b\ then the coarea formula holds: 

J? R (<p,x)dn u (x) = J dW(t) J dW-\u). (6.1) 

Proof. Obvserve that M = D U % U Z. Consequently, 
J J§ R ( 9 ,x)dn v (x) = J J§ R ( V ,x)dH»(x) 

M D 

+ J J§ R (<p,x)dn»(x) + J J§ R (tp,x)dH l/ (x). 

X z 
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Theorem 13.211 yields 

J J§ R ( v ,x)dU u (x) = J dU\t) j dH v - p (u). 

Theorem 14.11 implies 

j§ R (<p,x)dn u (x) = 0= f du\t) [ du v -\u) 



M V _1 (*)n3C 



Finally, since J7~ x ) = on Z, Theorem 15. II yields 

j J§ R {<p,x)dH"(x) = = J dU\t) J dU y -\u). 

z m (p-i(i)n0 

Combinig these relations, we obtain (16.11) . The theorem follows. □ 
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